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Abstract. It has been shown by Voiculescu and Biane that the analytic subordina- 
tion property holds for free additive and multiplicative convolutions. In this paper, we 
present an operatorial approach to subordination for free multiplicative convolutions. 
This study is based on the concepts of 'freeness with subordination', or 's-free inde- 
pendence', and 'orthogonal independence', introduced recently in the context of free 
additive convolutions. In particular, we introduce and study the associated multiplica- 
tive convolutions and construct related operators, called 'subordination operators' and 
'subordination branches'. Using orthogonal independence, we derive decompositions 
of subordination branches and related decompositions of s-free and free multiplicative 
convolutions. The operatorial methods lead to several new types of graph products, 
called 'loop products', associated with different notions of independence (monotone, 
boolean, orthogonal, s-free). We also prove that the enumeration of rooted 'alternat- 
ing double return walks' on the loop products of graphs and on the free product of 
graphs gives the moments of the corresponding multiplicative convolutions. 



1. Introduction 

Multiplication of free random variables X% and X2 with distributions Lt\ and 1x2, 
respectively, leads to the multiplicative convolution fj,i [x]£t2, introduced by Voiculescu 
[23] in the C*-algebra framework, which gives the distribution of the product of X\ 
and X 2 (the general case of measures with unbounded support was studied by Bercovici 
and Voiculescu [7]). 

Let .Mr + denote the set of probability measures on M + = [0, 00). If \i e A4r + , we 
can define 

(1.1) M*)=\ T^-dLi(t), zeC\R + , 

Jr+ l — zt 

which, in the case when fi has finite moments of all orders, becomes the moment gen- 
erating function ifj^(z) = ^™ =1 /i(X")2; n , where fi(X n ) are the moments of the unique 
functional \i : C[X] — > C defined by \i. In order to study //i[x]//2j Voiculescu introduced 
the S-transform of \i defined by S^{z) = (1 + z)ifj^ 1 (z)/z, where VC 1 ^) denotes the 
inverse of tpn( z ) with respect to composition. The key multiplicative formula for the 
S-transforms is given by S wi(1J (2) = S fll (z)S fl2 (z). 
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In our approach, a central role is played by the transform related to tp^z), namely 

where z e C\R + . Using this transform, Biane [8] has proved the subordination prop- 
erty for free multiplicative convolutions of probability measures on R + (and also for 
probability measures on the unit circle T). Earlier, the subordination property for free 
additive convolutions [22] was discovered by Voiculescu [24] for compactly supported 
measures on R, generalized by Biane [8] to arbitrary measures on R (see also [17] for a 
related approach). 

For instance, subordination for free multiplicative convolutions of probability mea- 
sures on R + says that for given fii,^2 e -Mr + \{c>o}, there exist analytic self- maps 771, r\ 2 
of C\R + , such that 

(!-3) Vi*mi»(z) = Vm(Vi(z)) = ffoM*))- 

One says that rj^^^ is subordinate to both r]^ and r)^ 2 , with rji and r/ 2 being the so- 
called subordination functions. These functions play a key role in the analytical study 
of free convolutions [3,4,10]. 

The functions rji and ^ 2 are unique and can be viewed as ^-transforms of certain 
probability measures on R + which are not concentrated at zero. This defines a binary 
operation on -Mr + \{c>o}, namely 

(!- 4 ) Vi(z) = ^ 20M1 (^), and rj 2 (z) = 77 W0W (2). 

The associated convolution /ii 0/i 2 , introduced in this paper, is called the s-free mul- 
tiplicative convolution and it plays the role of the multiplicative analog of the s-free 
additive convolution [i\ □ /j, 2 studied in [16]. 

The subordination formulas (1.3) are related to the so-called monotone multiplicative 
convolution of probability measures on R + , introduced and studied by Bercovici [6]. 
This convolution can be defined by the equation 

for z e C\R + , where /^i,/x 2 G -Mr + - Using this convolution and the s-free multiplicative 
convolution, we can write subordination equations (1.3) in the convolution form 

(1.6) |Ul@^2 =/*lO (Ai2 fJ-l) = (J-2 O (A*l A*2) 

for /i^/ia e -M M+ \{5 }. 

Further, we show that one can decompose both s-free and free multiplicative convo- 
lutions of compactly supported measures on R + which are not concentrated at zero in 
terms of simpler convolutions. For that purpose we introduce another convolution of 
/ii,/i 2 e -Mm + \{<5o} called the orthogonal multiplicative convolution, denoted ^iZ/i 2 . It 
plays the role of a multiplicative analog of the orthogonal additive convolution on R, 
introduced and studied in [16]. Using transforms, one can define this convolution by 

(1-7) W„W = 

"ill \ z ) 
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for z e C\R+. Using this convolution, we obtain a decomposition of the s-free multi- 
plicative convolution of \i\ \Z\ H2 in terms of an infinite sequence of alternating /ii,/J,2 if 
these are compactly supported. In view of (1.6), this leads to a decomposition of the 
free multiplicative convolution. 

In a similar way one can treat probability measures on the unit circle T, denoted 
Mj. In that case, rj^ is the integral over T and z lies inside the open unit disc D. One 
defines subordination functions for Hii^2 e -M.*, where M.* = A4j\{h : (J>(X) = 0}. 
Let us also observe that a number of results hold for distributions £ S, where 

S denotes the set functionals \i : C[X] — > C, which send 1 into 1. In that case, the 
functions 77^ are understood as formal power series. 

One of the main points of this paper is that our study of the relations between sub- 
ordination functions and their decompositions uses operatorial techniques. Namely, 
we construct bounded operators on Hilbert spaces, which correspond to all compactly 
supported convolutions which appear in the subordination equations and in the decom- 
positions of s-free convolutions. Our approach is based on the decomposition of the 
free product of Hilbert spaces (7Y,£) = (7ii,£i) * (7^2, £2), as (two different) orthogonal 
direct sums 

00 

(1.8) n = ® w (n_1) (0 ®n {n \i) 

n=l 

for each 1 e I = {1,2}, where 1 = 2 and 2 = 1, with liS n \i) denoting the subspace 
spanned by alternating tensor products of order n which do not end with a vector from 
7i G L : = 7Y t ©C£ t , following the original notation of Voiculescu [21] (see also [25]), and we 
set TL^\i) = C£ for each tel. By P L (n) we denote the canonical projection from Ti. 
onto H( n-1 )(i) ®H^ n \l) and we set to be the projection onto C£. Finally, we define 
the so-called vacuum state ip(-) = (■£,£) on B(TC). 

According to the above decomposition, we can represent bounded free random vari- 
ables in the form of 'orthogonal series' 

00 

(1-9) X b = Y,XXj), 

3 = 1 

where 1 e I and X t (j) = P L (j)X L P L (j). These series play a crucial role in our study of 
subordination for both free additive convolutions [16] and free multiplicative convolu- 
tions, studied in this paper. Let us remark that the 'orthogonal series' were introduced 
in a more general context of 'monotone closed *-algebras' of operators 'affiliated' with 
unital *-algebras [14]. In that approach, we used the 'monotone tensor product' (g) 
(reminding the von Neumann tensor product) to represent free random variables. 

Let H2 be (/^-distributions of bounded random variables X\ and X2, respectively, 
which are free with respect to ip. In the case of free additive convolutions, crucial was 
the decomposition of their sum in terms of 'additive subordination branches', namely 
Xi + X 2 = Bi + B 2 , with 



(1.10) 



B L = S° L dd + S? 
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where 

(1.11) Sr = J] X t (j) and Sr = J] X,(j), 

j odd j even 

for each tel. Although the subordination branch B L is a bounded operator on the 
'free Fock space' Ti for each i e I, it acts non-trivially only on its subspace, the 's-free 
Fock space' /C t , defined as 

(1.12) K L = 0^("- 1 )( i )©7YH(z), 

n odd 

and therefore, it is appropriate to view B b as an element of B(JC L ). 

For operatorial subordination in the mutliplicative case, we need to modify the op- 
erators S° dd , S c L vcn and take 

(1.13) R odd = S° L dd and Rl vcn = + S c L vcn . 

Moreover, in order to include all bounded positive random variables, we shall need the 
operation to be defined for all compactly supported /^i, /jl 2 e M.M.+ - For that purpose 
we set 

(1.14) n 5 = S^x) and S fi = S 

for compactly supported \i e A4r + , which turns out natural in our operatorial setting. 

Using these notations, and choosing, for simiplicity, an existence-type formulation, 
we can summarize our subordination result for positive operators as follows. 



Subordination for positive operators. If X 1 and X 2 are positive, then there 
exists a decomposition */X[X2'\fX[ = y/y Y^Jy, such that the <p- distributions of y and 
Y are ^ and fJ,2\Z\fJ>i, respectively, and the pair (y — 1, Y — 1) is monotone independent 
w.r.t. (p. Moreover, the if -distribution ofY agrees with that of i?J von R°2 A *J i?J vcn , and 
the pair (R^ — 1k 2 : -RT™ ~ l/c 2 ) * s s -free independent w.r.t. (<p, if)), where if) is the state 
associated with any unit vector £ e 

Moreover, we further decompose the 'positive subordination branches' 

(i.i5) v^r^V^r" 

using the notion of 'orthogonal independence'. These decompositions corresponds to 
decompositions of s-free convolutions in terms of orthogonal convolutions. In a similar 
way we obtain operatorial subordination results for unitary operators, or even more 
generally, for bounded operators. 

The main examples of operators related to subordination for the free additive con- 
volution, or, more generally, to various notions of independence X, are the adjacency 
matrices of subgraphs of the corresponding X-product graphs Q1IQ2 (in our study, X 
stands for orthogonal, comb, star, s-free or free). In fact, to each independence X, we 
can associate an additive and a multiplicative convolution, 



(1.16) 



Hi +1^2 and Hi x x jj 2 
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respectively. Recall that the additive X-convolution of spectral distributions of rooted 
graphs corresponds to the addition of 'monochromatic' X-independent adjacency matri- 
ces S L , t e I, and that, in turn, is related to the enumeration of rooted (i.e. root-to-root) 
walks on Q{LQ 2 (for details, see [2] and its references). 

In order to find a 'universal' multiplicative analog of this theorem, one needs to 
introduce a new concept of a product of Q\ and Q 2 for each X-independence, which we 
call the X loop product and denote Q{ItQ 2 . Roughly speaking, Q{IiQ 2 is obtained by 
adding colored loops to Q{IQ 2 in a suitable way. We assume that the product graph 
has a 'natural coloring', by which we mean that each copy of Q L is colored with color 
t, where cel. The procedure of adding loops is equivalent to replacing each S L by its 
'unitization' R L obtained from S L by adding some projection, in such a way that makes 
the pair (R\ — 1, R 2 — 1) X-independent. 

Quite naturally, in order to formulate our result, we shall use the formal power series 
rjz corresponding to the variable Z = R 2 R\, namely 



and interpret Nz(n) as the 'first return moments' in the state <p e associated with the 
root e. These moments are related to the enumeration of walks of the same class, which 
we find to be rooted alternating d-walks ('d-walk' is our abbreviation of 'double return 
walks originating with color 1') counted on different products. Finally, in view of the 
relation to independence X just mentioned, these moments agree with the 77-moments 
of the corresponding multiplicative convolutions /ii Xj /i 2 . It is not hard to see that 
this result can be generalized to the framework of random walks [26]. 

Multiplication theorem. Let Q\TtQ 2 be naturally colored and let A(Q{I e Q 2 ) = 
Ri + R 2 be the decomposition of its adjacency matrix induced by the coloring. Then 



where Z = R 2 R\ and D 2n (e) denotes the set of rooted alternating d-walks on Q\IiQ 2 of 
lenght In, where neN. Moreover, the pair (R x — 1, R 2 — 1) is X-independent. 

The paper is organized as follows. In Section 2, we introduce basic notions, including 
the s-free multiplicative convolution. In Section 3, we introduce and study the concepts 
of comb- and star loop products of graphs and find relations between rooted alternat- 
ing d-walks on these graphs and the monotone and boolean multiplicative convolutions, 
respectively. In Section 4, we introduce and study the orthogonal multiplicative con- 
volution. In Section 5 we define ans study the corresponding notion of the orthogonal 
loop product of rooted graphs. In Section 6 we show, by means of analytical methods, 
that the definition of the orthogonal multiplicative convolution can be extended to ar- 
bitrary measures on M + and T. The main operatorial results of the paper are contained 
in Sections 7 and 8, where we introduce and study operators on the free and s-free 
Fock spaces which are related to subordination for multiplicative free convolutions and 
their decompositions. Finally, in Section 9 we find a relation between free and s-free 



(1.17) 




(1.17) 
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multiplicative convolutions and the enumeration of rooted alternating d-walks on the 
free product of graphs and on the s-free loop product of graphs, respectively. 

Throughout the whole paper we understand that / = {1,2} and we adopt the nota- 
tion 1=2 and 2 = 1. 



2. Preliminaries 

By a non-commutative probability space we understand a pair (A,(p), where A is a 
unital algebra over C and ip is a linear functional ip : A — > C such that ip(l) = 1. If ^4 is 
a unital *-algebra and ip is positive (called a state), then (.4., y?) is called a ^-probability 
space. If, in addition, „4 is a C*-algebra, then (*4, tp) is called a C* -probability space. 
By the Gelfand-Naimark-Segal theorem, a C*-probability space can always be realized 
as a subalgebra of bounded operators on a Hilbert space 7i with a distinguished unit 
vector £, for which </?(a) = (a£, £) for a e A. 

By a random variable we will understand any element a of the considered algebra A. 
If „4 is equipped with an involution, then a random variable a will be called self-adjoint 
if a = a*. The ip- distribution of a random variable a is the functional /x a : C[X] — > C 
given by jtt a (l) = 1, /i a (X n ) = <p(a n ). In particular, if (A,ip) is a C*-probability space, 
then the distribution /i a of a self-adjoint random variable a e A extends to a compactly 
supported probability measure fj, on the real line. In that case we will often use the 
same notation \i for both the distribution of a and the associated compactly supported 
probability measure. 

Various notions of 'independence' X lead to several types of convolutions of distribu- 
tions (probability measures). If we have two random variables, Xi e A\ and X 2 e A2 
with (/^-distributions ji\ and fj, 2 , respectively, , where A± and A2 are X-independent sub- 
algebras of a noncommutative probability space (A,<p), then the additive convolution 
of /i and v associated with X-independence is the ^-distribution of the sum X\ + X 2 . 
In turn, the multiplicative convolution of /i and v is the ^-distribution of X 2 X\. In 
this paper we are interested in the free and s-free multiplicative convolutions associ- 
ated with free and s-free independence, respectively. In order to decompose them, we 
shall use the monotone multiplicative convolution introduced by Bercovici [6] and we 
will introduce the orthogonal multiplicative convolution associated with the notion of 
'orthogonal independence' [16]. 

Definition 2.1. Let fj, 1 , ji 2 e A^r + \{<5o} and let rj\ and r] 2 be the associated subordina- 
tion functions. The unique probability measures oi, a 2 e -Mr + \{5o} such that r\\ = r] a2 
and 7/2 = ?7 CTl will be called the s-free multiplicative convolution of [i\,[i 2 and [i 2 ,[i\, 
respectively. We set <J\ = \i\V\ A*2 and a 2 = fj, 2 In a similar way we define the 
s-free multiplicative convolution of fj,i, fj, 2 e M.*. 

The s-free multiplicative convolution defines a binary operation on both A1r + \{5o} 
and * • It can be seen that it is neither commutative nor associative. We will show 
later that it is related to the concept of 'freeness with subordination', or 's-free indepen- 
dence', introduced in [16] (a relation between freeness and monotone independence was 
also studied in [15]). Moreover, it allows us to rewrite (1.3) in terms of convolutions. 
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Proposition 2.1. If ^2 e A4r + \{5o}> ^en 

(2.1) |Ul@M2 = ^lO (a*2 0A*i) = ^2 O (a*i 0^2) 

is t/ie decomposition of the free multiplicative convolution corresponding to the subordi- 
nation equations (1.3). A similar result holds for jUi,^ e 

Proof. This fact is an immediate consequence of (1.3), in view of (1.5). ■ 

By a rooted set we understand a pair (X, e), where X is a countable set and e is a 
distinguished element of X called root. By a rooted graph we understand a pair (£?,e), 
where Q = (V, E) is a non-oriented graph with the set of vertices V = V(Q), and the set 
of edges E = E(Q) c: {(x, x') : x,x' e V} and e e V is a distinguished vertex called the 
root. We identify (x, x') with (x',x) since we consider non-oriented graphs, but when 
speaking of walks we find it convenient to say that (x, x') begins with x and terminates 
with x' . Another reason for using this notation for edges is that of main interest to us 
are graphs which have loops, i.e. edges of the form (x,x), where x e V. The notion 
of a rooted graph can be easily generalized to allow multiple edges. Formally, we then 
obtain a rooted multigraph, but we will still use the term 'rooted graph', or simply 
'graph' since all graphs will be considered to have a root. Moreover, we will very often 
omit the root in our notation and denote by Q the rooted graph (£/, e) if no confusion 
arises. Thus, in a graph Q we denote by n(x, x') the number of edges connecting x and 
x' (it may be zero). 

For (rooted) graphs we will also use the notation 

(2.2) V° = V\{e}. 

Two vertices x, x' e V are called adjacent if x and x' are connected by an edge, which 
we denote x ~ x'. The degree of x e V is defined by k(x) = ^ lx i^ x n(x,x'). A graph 
is called locally finite if k(x) < 00 for every x e V. It is called uniformly locally finite 
if sup{/t(x) : x e V} < 00. All graphs considered in this paper will be connected and 
uniformly locally finite. 

For x e V, let 5(x) be the indicator function of the one-element set {x}. Then 
{S(x), x e V) is an orthonormal basis of the Hilbert space 1 2 {V) of square integrable 
functions on the set V, with the usual inner product. The adjacency matrix A = A(Q) 
of a graph Q is a matrix defined by A(Q) XjX i = n(x, x') We identify A with the densely 
defined symmetric operator on l 2 (V) defined by 

(2.3) A{Q)5{x) = J] n(x,x')5(x') 

for x e V. Notice that the sum on the right-hand-side is finite since our graph is 
assumed to be locally finite. It is known that A{Q) is bounded iff Q is uniformly locally 
finite. 

By the spectral distribution of A(Q) in a state ip on h(V) we understand the measure 
\x for which 

(2.4) i>(A n ) = f x n /i(dx), n e N u {0}, 
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and by the spectral distribution of the rooted graph {Q, e) we understand the spectral 
distribution of A{Q) in the state <p e (-) = (-<H e ), <H e ))- 

A walk from v to v n on a graph Q is an alternating sequence of vertices and edges 
of the form 

(2.5) w = (v ,p 1 ,v 1 ,p 2 ,v 2 , ■ ■ . ,v n _i,(3 n ,v n ), 

with vertices vq,v\, . . . ,v n e V(Q) and edges Pi, p 2 , ■ ■ ■ , (3 n e E(Q), such that $ is an 
edge connecting i>j_i and V{. The lenght of w is given by l(w) = n. We allow t>j_i = 
in particular this happens if $ is a loop. 

A subwalk of w is a subsequence of w of the form 

(2.6) w' = (vi, Pi, v i+1 ,p i+1 , Pj-i, Vj), 

where ^ i < j ^ n. In the case when Q does not have multiple edges, we can identify 
w with the sequence (vq,Vi, . . . , v n ) with the understanding that i>j_i ~ Vi (in this case 
there is no confusion which edge connecting Vi_i and Uj is chosen) for 1 ^ % ^ n. 

The set of walks from v to v' will be denoted W(v,v') and we set W(v,v) = W(v). 
A walk from v = vq to v = v n , in which ^ # w for all 1 ^ i ^ n — 1 will be called 
an f-walk. Note that we do not require all vertices vq,v\, . . . , i> ra _i in an f-walk to be 
distinct. The set of f-walks from v to v will be denoted F(v). Subsets of W(v) and 
F(v) consisting of walks and f-walks of length n, respectively, will be denoted W n (v) 
and F n (v), respectively. 

A walk w can also be represented as a sequence of subwalks; in particular, an alter- 
nating sequence of edges and subwalks of the form 

(2.7) W = (W , P 1 ,W 1 , p 2 , W 2 ,..., W r -i, Pr, W r ), 

where w is a walk from v to some v^, Pi is an edge connecting with W\ 
is a walk from fj(i)+i to some v^, P 2 is an edge connecting v^ 2 ) with i>j(2)+i, and so on, 
finally, w r is a walk from i>i( r )+i to Uj( r +i) = v n . We shall also use similar representations 
of walks which begin or end with an edge. Representing a walk in terms of edges and 
subwalks is very convenient when the subwalks W\ , w 2 , . . . , w r are of special type (for 
instance, are f-walks). 

A graph Q whose edges are colored with colors from the set / will be called an I- 
edge-colored graph. In our case, we will always assume that / = {1,2}, so we can also 
use the term '2-edge-colored graph'. A walk in such a graph will be called alternating 
if the colors of its edges alternate. 

Definition 2.2. Let Q be /-colored. A walk w e W(v) on Q will be called an alternating 
double return walk, or simply an alternating d-walk, if it is alternating, begins with 
an edge of color 1, and can be represented as a pair (ui,u 2 ) of subwalks, such that 
u\,u 2 e F(v). The set of alternating d-walks from v to v will be denoted D(v ), and we 
set D m (v) = W m (v) nD(v). 

Edge-coloring by a two-element set is natural in the case of many products of rooted 
graphs. In particular, the products of rooted graphs of type Q{IQ 2 have the property 
that every edge of the product graph belongs either to a copy of Q\, or to a copy of Q 2 . 
Thus we can color the edges of all copies of Q L with color t, where t e I, in which case 
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we will say that Q{IQ 2 is naturally colored. In turn, we will say that the product Q\T f Q 2 
is naturally colored if its coloring is inherited from Q{IQ 2 . Let us also remark that our 
condition that an alternating walk should begin with an edge of color 1 (and thus end 
with an adge of color 2) is caused by the fact that we want to identify the given walk 
w and its 'reverse' obtained by reversing the order in (2.12). 

We end this Section with elementary formulas for the 'first return moments' of a ran- 
dom variable X, by which we understand the coefficients of the formal power series r] X 
associated with the distribution of X. If X is a random variable in a non-commutative 
probability space (A,tp), then the 'moment generating function' and the 'first return 
moment generating function', respectively, associated with the (/^-distribution fi of X, 
are given by formal power series 

oo , / \ oo 

(2.8) ^) = J>pr)(nK\ Vx (z) = 1 ^ Z) =^N x (n)^ 

n=l tX\ J n=l 

where the numbers Nx(n), n eN, will sometimes be called 'first return moments of X\ 
Below we give a convenient algebraic formula for these 'first return moments'. For 
that purpose, let us take the extension of (A, (p) to a larger noncommutative probability 
space, namely, the free product with identified units A * C[P], where P 2 = P, with the 
state given by the linear extension of 

(2.9) <p(P) = 1, <p( Wl Pw 2 ) = <p( Wl )<p(w 2 ) 

for any w\,w 2 e A* C[P], where, slightly abusing notation, we also denote the new 
state by ip [13]. In particular, in the C*-algebra context, P can be interpreted as the 
projection onto C£, where £ is the cyclic unit vector of the GNS triple associated with 
(A,<p). 

Proposition 2.2. Let X be a random variable in a noncommutative probability space 
(A, if). Then 

(2.10) N x (n) = i P (X(P L X) n - 1 ) 
for n eN, where P L = 1 — P. 

Proof. This is a straightforward consequence of (1.2). ■ 

Finally, let us prove an elementary fact about a relation between the 'first return 
moments' of products of 'monochromatic' adjacency matrices in certain /-colored graphs 
and the cardinalities of the sets D 2n (e). 

Proposition 2.3. Let (Q, e) be an I -colored uniformly locally finite graph, and let 
A(Q) = A\ + A 2 be the decomposition of A(Q) induced by the coloring. If the set 
of rooted alternating f-walks of even lenghts is empty, then 

(2.11) N z (n) = \D 2n (e)\, 

where Z = A 2 A\ and the numbers Nz(n) correspond to the state <p e (-) = (■^( e )^( e ))- 

Proof. If we take Z = A 2 A X in Proposition 2.2, then N z (n) is equal to the number of 
rooted alternating walks of lenght 2n, which begin with an edge of color 1 and such that 
intermediate returns to the root can occur only after odd numbers of steps. If there 
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are no rooted alternating f- walks of even lenght, then the first return to the root must 
occur after an odd number of steps and thus the second return to the root must occur 
after an even number of steps. Therefore, if we consider rooted alternating walks which 
have an even number of steps and such that intermediate returns occur only after an 
odd number of steps, these are d-walks. This proves our assertion. ■ 

Let us note that all graph products considered in this paper satisfy the assumptions 
of Proposition 2.3, and therefore, walk-counting is always reduced to rooted alternating 
d-walks. 



3. Comb and star loop products of graphs 

We begin with recalling the notions of the additive [18] and multiplicative [6] mono- 
tone convolutions and show that the 77-moments of the latter are related to the enu- 
meration of alternating d-walks on a new version of the comb product of graphs called 
the 'comb loop product'. 

Definition 3.1. Two subalgebras A±,A2 of a unital algebra A are monotone indepen- 
dent with respect to a normalized linear functional ip on A if 

<p(wiaib) = (p(b)(p(wiai) and (p(ba 2 w 2 ) = (p(b)p(a 2 w 2 ), 
(p(wiaiba 2 w 2 ) = (p(b)(p(wiaia 2 w 2 ), 

whenever cii,a 2 e A±, b e A 2 and w±,w 2 are arbitrary elements of the unital algebra 
alg(^4i,^4 2 ) generated by Ai and A 2 . In particular, we will say that the pair (a,b) 
of elements of A is monotone independent w.r.t. ip if the algebras generated by these 
elements are monotone independent. 

Suppose (^-distributions of random variables a± and a 2 are fj,i and fj, 2 , respectively. If 
the pair (a l5 a 2 ) is monotone independent w.r.t. ip, then the ^-distribution of the sum 
01 +a 2 is the monotone additive convolution jii t> /j, 2 . In turn, if the pair (a x — 1, a 2 — 1) 
is monotone independent w.r.t. ip, then the (^-distribution of the product a 2 (i\ is the 
monotone multiplicative convolution [i\ Q /i 2 . The corresponding formal power series 
satisfy (1.5). If fj,i, /i 2 e A4r + or fi\,fi 2 e A4j, then (1.5) is understood in terms of 
analytic self-maps of C\M + or T, respectively. For details, see [19] (additive case) and 
[6] (multiplicative case). 

It is known that the monotone independence can be associated with the comb product 
of graphs [1]. Let us recall the definition of the comb product and follow up with the 
corresponding loop product. 

Definition 3.2. The comb product of rooted graphs (Qi,ei) and (Q 2 ,e 2 ) is the rooted 
graph (Qi t> Q 2 , e) obtained by attaching a copy of Q 2 by its root e 2 to each vertex of 
Qi, where we denote by e the vertex obtained by identifying e\ and e 2 . If no confusion 
arises, we denote the comb product by Q\ o Q 2 . If we identify its set of vertices with 
Vi x V 2 , then its root is identified with e\ x e 2 . 
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Definition 3.3. Suppose the edges of the comb product (Q 1 t> Q 2 ,e) are naturally 
colored. The comb loop product of rooted graphs (Qi,ei) and (Q 2 ,e 2 ) is the rooted 
graph (Qi \>- e Q 2 ,e) obtained from (Q\ t> Q 2 ,e) by attaching a loop of color 1 to each 
vertex but the root of each copy of Q 2 . 

Let us justify the above definition. We follow the observation made by Bercovici 
[6] that in order to introduce a multiplicative convolution associated with monotone 
independence, one needs to 'unitize' the usual monotone independent variables. It has 
been shown in [5] that one of the possible choices is to take variables of the form 

(3.1) i?i = ai <g> P 2 + li ® P 2 and R 2 = l 1 ®a 2 , 

where a L , 1 L e B(7i L ), with (Ti t , £ t ), being Hilbert spaces with distinguished unit vectors 
and identity operators l t , and tel. Then the moments of the product R 2 Ri in the 
state associated with the vector £1 (x) £ 2 agree with the moments of /xi Q /i 2 . 

If a t is taken to be the adjacency matrix of a (uniformaly locally finite) graph Q L , 
1 e I, then the term I4 (x) corresponds to glueing a loop of color 1 to each vertex 
but the root of each copy of Q 2 . Thus, R\ is obtained from the usual comb product 
adjacency matrix S\ by adding a projection L±, wheras R 2 = S 2 , where 

(3.2) S ± = ai(g)P 2 , L x = l 1 (g)P 2 L 

and Ri, R 2 are the adjacency matrices of 'monochromatic', 1 and 2-colored subgraphs, 
respectively. This realization can be generalized to arbitrary random variables in non- 
commutative probability spaces, using the extensions of functionals given by (2.9). 

Remark 3.1. It is not difficult to show that the comb loop product of graphs Q\ and Q 2 
is isomorphic to the usual comb product of larger graphs Q\ and Q 2 . However, in order 
to keep the same coloring of the product graph (wich is needed for counting alternating 
walks) one has to use two colors for Q 2 . Namely, Q 2 is obtained from Q 2 by glueing a loop 
of color 1 to each vertex but the root of Q 2 . Then we indeed have Q\ t> e Q 2 = Q 1 o Q 2 , 
but usefulness of this relation seems to be limited (roughly speaking, we have a simpler 
product but a more complicated coloring). 

Example 3.1. Let us consider the example of the comb loop product of graphs given 
in Figure 1. We choose both Q\ and Q 2 to have a loop at the root (to distinguish the 
loops, we draw loops of color 1 smaller than loops of color 2, which is helpful in the 
enumeration of alternating walks). Then Qi t>^ Q 2 has loops of two types: loops whose 
origin can be traced back to graphs Q\ and Q 2 (these are all loops which are at the 
glueing points), or are added to the usual comb product in the process of forming the 
comb loop product. Denote by (v, v\ the loop from v to v of color 1. The following 
walk is a rooted alternating d-walk: 

(/3i, wx, fa, w 2 , As, u> 3 , /3 4 , Wi) e D 8 (e), 

where the sequence of edges of color 1, 

(Pi, fo, 03, 0a) = ((e, x), (x, x'), (x 1 , x), (x, e)) 

forms a rooted f-walk of color 1, which is interlaced with the sequence (wi, w 2 , W3, W4) 
of loops of color 2: w\ = w 3 = (x, x) 2 e Fi(x), w 2 = (x 1 , x') 2 e Fi(x') and w 4 = (e, e) 2 e 
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Figure 1 . An example of the comb loop product 

F 1 (e). Note that the first return to e occurs at the end of fa, whereas the second return 
occurs as the end of w±. When any of the loops, w i: 1 ^ % ^ 4, is replaced by an 
arbitrary alternating f-walk from the corresponding vertex to itself which begins and 
ends with an edge of color 2, we still get a rooted alternating d-walk. The simplest 
examples of a rooted alternating d-walk are of course given by: (fa, Wq) e .02(e), where 
fa = (e, e)i, w = (e, e) 2 and (fa, w ± , fa, w 4 ) e D 4 (e). 

Theorem 3.1. The Multiplication Theorem holds for monotone independence, the as- 
sociated loop product Qi o- e Q 2 , and the multiplicative convolution [i\ Q /i 2 . 

Proof. From Definition 3.3 and (3.1) it follows that 

z = A(g 1 1>, g 2 ) = r 1 + r 2 , 

where the pair (R 1 — 1,R 2 — 1) is monotone independent w.r.t. ip e , the state associated 
with e and R L , 1 e /, is the adjacency matrix of the ^-colored subgraph. Moroever, 
by Proposition 2.2, Nz(n) is the number of alternating d-walks on Qi t>£ Q 2 of lenght 
2n. This proves the first equation in (1.13) for monotone independence. Now, from 
equation (1.5), we obtain the combinatorial formula 

n 

N^ 2 (n) = J] iV M1 (r) J] N IJi2 (k 1 )N IX2 (k 2 ) . ..N, 2 (k r ), 

r—l ki+k2 + ...+k r =n 

where it is tacitly assumed that the indices k±,k 2 , . . . ,k r are positive integers and num- 
bers A^(n) are coefficients of rj^ treated as formal power series 

00 

= J] N,(n)z n . 

n=l 

The above combinatorial formula allows us to find a correspondence between these 77- 
moments and rooted alternating d-walks on Q\ x> t Q 2 . Let us observe that (r) is the 
number of rooted f-walks of lenght r on Q 1 and N^ 2 (k) is the number of rooted f- walks 
of lenght k on Q 2 . On the other hand, recall that in the comb loop product Q\ c>£ Q 2 
there is only one copy of Q\ (with e\ identified with the root e of the product graph), 
with a copy of Q 2 attached by its root to every vertex x of Q\ (the vertex x becomes then 
the only common vertex of this copy of Q 2 and the original copy of Qi). In addition, 
loops of color 1 are glued to all vertices but the roots of these copies of Q 2 . Therefore, 
each rooted alternating d-walk on Q 1 Q 2 consists of a sequence of edges of Q\ which 
forms an f-walk of color 1, namely c = {vq,V\, v 2 , . . . ,v r -i,v r ) e F(e), and alternating 
f-walks q e F(vi), i = 1, . . . , r, with the first edge of color 2, attached to every vertex 
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of c. Note that the number of alternating f-walks from Vi to Vi with the first edge of 
color 2 is equal to the number of walks from e 2 to e 2 in graph Q2 - the only difference 
is that if the latter has lenght ki, the fomer has lenght 2/cj — 1 since each edge must be 
followed by a loop of color 2. Thus, the contribution from each product of the form 

N^(r)N^(h)N^(k 2 )...N^(k r ) 

to the RHS of the formula for iV Ml Q M2 (n) is equal to the number of all such f-walks c of 
color 1 which have r edges, interlaced with r alternating f-walks of lenghts 2ki — 1, 2k 2 — 
1, . . . , 2k r — 1. The summation over k\ + k 2 + . . . + k r = n indicates that all alternating 
f-walks involved must have 2n — r edges together. The summation over 1 ^ r ^ n gives 
exactly the cardinality of D 2n (e), which finishes the proof. ■ 

Let us observe here that (3.1) can be generalized to arbitrary random variables in 
noncommutative probability spaces with distributions /xi,/x 2 e E if one treats Pi and 
P 2 as idempotents and uses extensions of states ipi, ip 2 given by (2.9). Then, the com- 
binatorial formula in the proof of Theorem 3.1 remains valid for \i\,\i 2 6 E. 

Example 3.2. If /ii,/i 2 e E, the lowest order moments of iV W Q M2 are given by 

N^ 2 (l) = Ar m (l)AUl), 

JWP) = W W (2X(1) +^(1)^(2), 

JW(3) = ^(3)^(1) +2^(2)^(2)^(1) + ^(1)^(3), 
AW(4) = iV M1 (4X(l)+37V /ii (3)^ 2 (2)7V2 2 (i ) + 7V w (2)7V2 2 ( 2 ) 
+2JV W (2)A^ 2 (3) JV M (1) + iV m (1) JV M (4). 

Let us apply these formulas to the graph Q 1 c=- e Q 2 in Fig.l. Let [i\ and fi 2 be the spectral 
distributions associated with Q\ and Q 2 , respectively. By counting f-walks on Q\ and 
G 2 , we easily get N^l) = iV w (2) = iV w (4) = 1, 7V m (3) = and N^(l) = N^ 2 (2) = 1, 
A^ 2 (3) = 0, A^ 2 (4) = 2. In turn, using the above formulas and Theorem 3.1, we get 
D 2 (e) = 1, Di{e) = D 6 (e) = 2 and D 8 (e) = 4, which can also be verified directly by 
counting comb walks on Q\ t>£ Q 2 . 

For completeness, let us also briefly discuss the case of the star product of graphs 
Q\ * Q2 and find its relation to the boolean mutliplicative convolution. 

Definition 3.4. The star product of rooted graphs (£/i,ei) and (Q 2 ,e 2 ) is the rooted 
graph (Q\ *{? 2 , e) obtained by glueing Q\ and Q2 at their roots and taking this common 
root to be the root of the product. The star loop product (Q 1 * £ Q 2j e) is obtained from 
naturally colored (Qi * Q 2 , e) by adding loops of color 1 to every vertex but the root of 
where tel. 



Recall here the definition of the boolean multiplicative convolution of distributions 
H\iH 2 £ S [5]. Treating rj^z) as a formal power series, we define 

(3.3) p,(z) = 



14 



R. LENCZEWSKI 




Figure 2. An example of the star loop product 

which allows us to define the boolean multiplicative convolution of pi,p 2 e £ by the 
formula 

(3-4) p mSM2 (z) = p m (z) Pll2 0) . 

This formula defines H as a binary operation on Mj [11]. However, if pi,p 2 e -Mr + , 
then the RHS does not always give the p-transform of a probability measure on R + [5] . 
Nevertheless, the well-known relation [8,24] 

(3-5) PmEw(z) = Pi( z )P2(z) 

where p L (z) = r} L (z)/z, t e I, can still be written in terms of convolutions as 

(3.6) pi \X\(J> 2 = (pi 0P2) (P2 0Pi)- 

which was our motivation to use the symbol for the s-free convolution (an analogous 
formula holds for pi, p 2 e .M*). 

If a L is the adjacency matrix of a graph £/ t , where 1 e I, then, using a similar reasoning 
as in the case of the comb loop product, we introduce operators of the form 

i?i = ai ®P 2 + l®P^ 

R 2 = P l ®a 2 + Pj L ® h, 

which can be viewed as the adjacency matrices of the 1- and 2-colored subgraphs of the 
product graph (Qi *£ Qi,e). Moreover, these are exactly the operators which give the 
moments of the boolean multiplicative convolution of spectral distributions of Q\ and 
G2 since the pair — 1, R 2 — 1) is boolean independent w.r.t. <p e . An example of a 
star loop product is given in Fig. 2. 

Let us also note that formulas (3.3)-(3.4) can be used if ai,a 2 are random variables 
in arbitrary noncommutative probability spaces - this is obtained if one uses extended 
states (2.9). This gives a realization of the boolean multiplicative convolution of arbi- 
trary distributions pi,p 2 e E. 

Proposition 3.1. The Multiplication Theorem holds for boolean independence, the as- 
sociated loop product Qi *£ Q 2 , and the multiplicative convolution pi \*\ p 2 . 

Proof. The proof is similar to that of Theorem 3.1 and is based on the easy combi- 
natorial formula 

j+k=n+l 

where the sum runs over positive indices j, k, obtained from (1.2) and (3.4). Details 
are left to the reader. ■ 
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Example 3.3. If jUi,/x 2 e E, the lowest order 'first return moments' A^ lg]M2 (n) are 
given by 

= AWl)AU(l), 
N, im2 (2) = ^(1)^(2) + ^(2)^(1), 

N, im ,M = N^(1)N^(3) + AT m (2)AW2) + AW3)AU(l), 

^ lS ^ 2 (4) = iV m (I)Ar M2 (4) + Ar Ml (2)Ar M2 (3) + Ar Ml (3)Ar M2 (2) + Ar Ml (4)Ar M2 (I). 

These formulas can be used to count rooted alternating d-walks on the product graph 
given in Fig. 2. Using Proposition 3.2 and the values of N^^j) and N fl2 (k), given in 
Example 3.1, we get D 2 (e) = D A (e) = D e (e) = 1 and D 8 (e) = 3, which can be also 
verified by direct computations. 



4. Orthogonal multiplicative convolution 

For given positive bounded random variables Xi and X 2 with distributions fj,i and fi 2) 
we want to find a positive random variable Z on some Hilbert space whose distribution 
is given by the multiplicative analog of the orthogonal convolution, denoted /xiZ/u 2 . 

Let us first recall the notion of orthogonal independence [16]. 

Definition 4.1. Let (A,<p,ip) be a unital algebra with a pair of linear normalized 
functionals and let A\ and A 2 be non-unital subalgebras of A. We say that A2 is 
orthogonal to A± with respect to (ip, if)) if 

(i) (f(bw 2 ) = ip{wib) = 0, 

(ii) (f(w 1 a 1 ba 2 W2) = tp(b) ((p(wia 1 a 2 w 2 ) - (p(w 1 a 1 )ip(a 2 w 2 )), 

for any a\,a 2 e Ai, b e A 2 , and any elements wi,w 2 of the unital algebra alg(„4i,»4. 2 ) 
generated by A\ and A 2 . We say that the pair (a, b) of elements of A is orthogonal 
with respect to (ip, ip) if the algebra generated by b e A is orthogonal to the algebra 
generated by a e A. 

Definition 4.2. Let (^4, ip, ip) be a unital algebra with a pair of normalized linear 
functionals and let Ai,A 2 be a pair of (in general, non-unital) subalgebras of A, such 
that A 2 is orthogonal to A\ w.r.t. (<p,ip). Let ai,a 2 e A be random variables with 
(^-distribution ji\ and -^-distribution fi 2 , respectively. The orthogonal multiplicative 
convolutionof Hi and fi 2 , denoted //iZ/i 2 , is the distribution of a 2 ai, where (a\ — 1, a 2 — 1) 
is orthogonal w.r.t. V0- 

In order to find a Hilbert space realization of the orthogonal multiplicative convolu- 
tion of compactly supported probability measures, recall the definition of the orthogo- 
nal product of two Hilbert spaces with distinguished unit vectors (Hi,£i) and (7i 2 ,£ 2 ). 
Thus, the orthogonal product of (Hi,£i) and (7Y 2 ,£ 2 ) is the pair (7i,£), where 

(4.1) H = C£®H 1 ®(H 2 ®H° 1 ), 

with = Hi ©C£j and where £ is a unit vector. We denote it by (Hi,£i) \- (7i 2 ,£ 2 ), 
or simply Hi h H 2 , if no confusion arises. It has been shown in [16] that H gives a 
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Hilbert space realization of orthogonal random variables. For that purpose we defined 
an isometry U : H — > Hi ® H 2 given by 

(4.2) C/(0 =6® 6, ^i) = /^i®6, U(h 2 ®h 1 ) = h 1 ®h 2 

for any /ii e 7^5 and /12 e 7^2- Using U, we define faithful non-unital *-homomorphisms 
r t : > B(H) by 

(4.3) r^aO = U*( ai ®P 2 )U, r 2 (a 2 ) = U* (P, 1 ® a 2 )U, 

where P±, P 2 are the projections onto C£i and C£ 2 , respectively. Then the pair 
(ti(cii) , T 2 (a 2 )) is orthogonal w.r.t. (tp,ip), where ip and if) are states on B(H) asso- 
ciated with the unit vector £ and any unit vector £ e respectively. For details, see 
[16]. 

For simplicity, we will identify vectors from H with their images in Hi ® H 2 under 
the isometry U as well as bounded operators a, b with Ti(ai), r 2 (a 2 ), respectively, and 
we will carry out computations in Hi ®H 2 . 

Proposition 4.1. Let a x e B{H\) have tpi- distribution and a 2 e B(H 2 ) - <p 2 - 
distribution \x 2 . Then there exist a Hilbert space H with unit vectors £ and (, and ran- 
dom variables R\,R 2 e B(H) such that R\ has ip- distribution fii, R 2 has if)- distribution 
fj, 2 and R 2 R\ has <p- distribution [i\/-\i 2 , where ip and if) are states associated with £ and 
(, respectively. 

Proof. Let H = Hi \- H 2 with the canonical unit vector £ = d®^ 2 and let ( = (i®£, 2 , 
where £i is an arbitrary unit vector from H\. Define 

i?i = a x ®P 2 + h®P 2 , 
R2 = Pi®a 2 + P 1 ®\ 2 , 

where P^ = 1; — Pj and Pj is the canonical projection onto and lj denotes the 
identity in B(Hi), i = 1,2. Therefore, denoting 1 = \± ® 1 2 , we get 

Ri-l = (a 1 -h)®P 2 , 
R 2 -l = P 1 1 ®(a 2 -1 2 ), 

which, due to the appropriate tensor product form, implies that the pair (Ri — 1,R 2 — 1) 
is orthogonal w.r.t. (<p,ip), where tp is the state associated with £ and if) is the state 
associated with (. Moreover, elementary calculations show that the 92-distribution of 
Ri is iii and that the ^-distribution of R 2 is /i 2 , which completes the proof. ■ 

Corollary 4.1. Let fi\, fi 2 be compactly supported probability measures on R+. Then 
there exist bounded positive random variables R\, R 2 on some Hilbert space H and unit 
vectors £, £ e H, such that Ri has p- distribution Hi, R 2 has ip- distribution // 2 and 
\[P 2 R\\[P 2 has tp- distribution /iiZ/i2- 

Proof. There exist positive operators Xi e B(L 2 (R + , fii)) and x 2 e B(L 2 (R + , /i 2 )) 
(standard multiplication operators) which have distributions /ii and /i 2 , respectively, 
with respect to the states ipi and p> 2 associated with the constant functions equal to 
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one, hi and h 2 , respectively. We can then use the framework of Proposition 4.1 to get 
self-adjoint bounded random variables on H = L 2 (IR, /ii) (x) L 2 (R, /i 2 ) of the form 

R[ = (x 1 -l 1 )®P 2 , 

R' 2 = P 1 ± ®(x 2 -1 2 ), 

which are orthogonal w.r.t. (<f, ip), where if is associated with vector £ = hi (x) h 2 and ip 
is associated with ( = f<S>h 2 , where / is any function such that $ R f(x)d/ii(x) = and 
$ R f 2 (x)dfi 1 (x) = 1. Now, let R± = R[ + 1 and R 2 = R' 2 + 1. It follows from the proof of 
Proposition 4.1 that R± and R 2 have distributions yUi and fi 2 , respectively. Moreover, 
both are positive operators and thus, in particular, y^R^ exists. In fact, 

and ^JR7 2 R\\/R~ 2 is a bounded positive operator on H. Moreover, ^JR 2 Rl^fR 2 has the 
same (/^-distribution as R\R 2 , which is obtained from the following calculation: 

v?((V^iV^r) = vKV^i^r^iV^) 

= (fi ® <P2)((VP~i ® h^RiR^RiiPi <g> 1 2 )) 

= ^((iW 1 ). 

By using the adjoints, we can interchange the order of R\ and i? 2 under the symbol 
if. Therefore, the probability distribution of ^/R^Ri^/R^ extends to the compactly 
supported measure /iiZ/x 2 on R + . ■ 

We will now compute the (^-distribution of the product Z := R 2 R±, where R\ and 
R 2 are given by Proposition 4.1, and derive a formula for r\ z . In this context, let us 
observe that one can generalize Proposition 4.1 to arbitrary random variables in non- 
commutative probability spaces since formulas for Ri and R 2 remain valid provided one 
takes extensions of states given by (2.9) and uses Proposition 2.2. Therefore, further 
developments in this Section will hold for arbitrary random variables in noncommu- 
tative probability spaces. This corresponds to taking distributions n e S and the 
corresponding formal power series rj^ and p^. 

In order to compute the ^-distribution of Z, let us decompose Z as 

Z = Z\ + Z 2 + Z% + .£74, 

where 

Z 1 = P^ai ® a 2 P 2 , Z 2 = ® a 2 P 2 L , Z 3 = P^ (x) P 2 , Z 4 = P,® P 2 L . 

In order to compute the 'first-return moments' Nz(n) in the state f, we need to compute 
the mixed (^-moments of Z^s. However, it turns out that there are only two types of 
them which give a non-vanishing contribution. We compute them in Propositions 4.3- 
4.4. 

Proposition 4.2. For r odd and j\, j 2 , . . . , j r ^ 1, we have 

f(Z :i /r/^Zr • . . Z{ r ) = N^N^h) . . . N^k^), 
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where m = 1 + ^ kodd jk o,nd all hi 's are equal to 1, except 

kjl+l = h + 1) kjl+ji + l = J4 + lj • • • ! ^jl+i3 + ---+ir-2 + l = Jr-1 + 1) 

and where the 'first return moments' N aL (s) refer to the state ip L , tel. 
Proof. We have 

tpiZsZfZi'Zi 3 . . . Z?) = ^{P^iPta^PtiPt^r ■ ■ ■ PiiP^aiYn 

V2{P2{a2P2) n {a 2 P^(a 2 P 2 y*- 1 . . . (a 2 P 2 )>) 

= ¥ »i(ai(P i L aiy' 1+ « + - + >) 

V2{{a2P2) n {a 2 P 2 L y*{a 2 P 2 y* . . . (a 2 P 2 )>) 

= N ai (l+j 1+ j 3 + . . . + j r )(JV« a (l)) il+W '- 1)+ - +W, - 1) 
iV a2 ( j2 + 1) iV a2 (j 4 + 1) . . . AT a2 (j"r-l + 1), 

using the following properties: 

P101P1 = Vi(ai)Pi = ^(1)^1, P20-2P2 = ^ 2 (a 2 )P 2 = iV a2 (l)P 2 . 
Setting £4 = /c 2 = . . . = fcjj = 1, kj 1+ i = j 2 + 1, kj 1+2 = . . . = kj 1+ j 3 = 1, kj 1+ j s = 1, 

^Jl+i3 + l = J4 + I5 • • ■) %i+i3 + ---+ir-2 + l = Jr-1 + 1) %i+i3 + ---+Jr-2+2 = • • • = = 1, We 

get the desired result. ■ 

Proposition 4.3. For r even and j±,j 2 , . . . , j r ^ 1, we /iave 

^(Z 3 Zf Z? Z? . . . Z( r ) = ^(m)^^) . . . AU^-i), 
where m = 1 + Sfceven Jfc on< ^ a ^ ' s are e ? ua ^ ^° 1> except 

kj 2 + l = J3 + 1) kj 2+ j i+ \ = J5 + 1, . . . , kj 2+ j i + ,,, + j r _ 2 + i = j r -i + 1, 

and where the 'first return moments' iV 0t (s) refer to £/ie stote </? t; te/. 
Proof. Using similar arguments as in the proof of Proposition 4.3, we get 

zf/f . . . z?) = ^{p^p^p^rp^p^r ■ ■ ■ (AVn 

MP2(a2P 2 ± ) n (a2P2y 2 (a2P 2 1 y 3 . . ■ (a 2 p 2 n 

= (^ 1 (a 1 (P 1 ± a 1 ) J2+ ^ + - + >) 

^((a 2 P 2 ^y i (a2P2) j2 (a 2 P 2 ± ) j3 . . . (a 2 P 2 )>) 

= AT ai (l + j 2 + J 4 + . . . + j r )^a a (l)) ( * 2 - 1) + C * 1 - 1) + - + tfp - 1) 
N a2 ( 3l + l)iV a2 (j 3 + 1) . . . N a2 (jr-1 + 1). 

Setting h = ji + 1, k 2 = . . . = k j2 = 1, k j2+1 = j 3 + 1, & j2+2 = . . . = k j2+j4 = 1, 

kjl+ji + l = J5 + 1) • • •! %2+J4 + ...+ir-2 + l = Jr-1 + 1) %2+J4 + ...+Jr-2+2 = • • • = ^m-l = 1, We 

get the desired result. ■ 

Theorem 4.1. If the tp -distribution of a 2 is not concentrated at zero, then the formal 
power series corresponding to random variables ai,a 2 and Z satisfy the relation 

(4.4) „ z (z) = ^Ml. 

Va 2 0) 
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If the distribution of a 2 is concentrated at zero, then rjz(z) = zfi 1 (X), where fii is the 
(p -distribution of a±. 

Proof. Suppose the (^-distribution of a 2 is not the Dirac delta. Then, writing the 
RHS of equation (4.4) as D(z) = Yin^i^z 11 , we need to show that d n = r)z{n) for 
n eN. We clearly have 

N Z {1) = <p(Z 3 ) = v^i(ai) = N ai (l) = d l . 
Now, for n > 1, it holds that 

n 

d n =Y, N -M) 2 N ai (k 1 )N a2 (k 2 )...N a2 (k m _ 1 ), 

m=l fci+fc2 + ...+fc m _i=n— 1 

where all k^s are assumed to be positive integers, Using Proposition 2.2 and the explicit 
form of Zj's, we get 

N z (n) = ^{Z{P^Z) n - 1 ) 

= cp(Z 3 (Z 1 + Z 2 ) n - 2 Z 1 ) 

= S S 2 <p(z 3 z»z»...z»;z 1 ), 

where ji,j 2 , ■ ■ ■ ,j r are positive integers. Note that on the RHS of the formula for d n 
we have the sum of products 

N ai (m)N a2 (ki)N a2 (k 2 ) . . . N^k^), 

where 1 ^ m ^ n and (k±,k 2 , . . . , k m _i) is a sequence of positive integers which add 
up to n — 1. Now, every product obtained in Propositions 4.3-4.4 is of exactly this 
form, with indices (k±, k 2 , . . . , k m _ 1 ) satisfying the conditions just mentioned. Moreover, 
all these products are different from each other and exhaust all possibile values of 
(m, ki, . . . , k m _i). In fact, there are two disjoint cases distinguished: Proposition 4.3 
covers the case ki = 1 and Proposition 4.4 - the case ki 1. In the first case, even 
ji's determine the subsequence of (k 1: k 2 , . . . , & m _i) consisting of numbers which are 
greater than 1 and odd j'j's determine their positions in the sequence by determining 
all fcj's which are equal to 1. In the second case, odd jV s determine the subsequence of 
(ki, k 2 , . . . , fc m _i) consisting of numbers which are greater than 1 and even j^s determine 
their positions in the sequence by determining all /cj's which are equal to 1. Therefore, 
both cases give together every tuple (m, k±,k 2 , . . . , k m _i) under consideration exactly 
once. This proves (4.4) in the case when the distribution of a 2 is not concentrated at 
zero. Finally, if the 92-distribution of a 2 is So, then ^2(^2) = S n ,o, which easily gives 
Nz(n) = N ai (l)5 n: i, where 5ij is the Kronecker symbol, and since iV ai (l) = ipi(a\) = 
fj,i(X), we have rjz(z) = zfj,i(X). This completes the proof. ■ 

This shows that Eq.(1.7) holds for the orthogonal multiplicative convolution of dis- 
tributions yUi,yU 2 e provided /i 2 ¥= S , where the functions r\ are treated as formal 
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power series. In turn, explicit computations of distributions based on Definition 

4.2 give 

(4.5) 5qA^i = 5 and /xZ5 = S^ X ) 

for any /ieS. Let us also remark that Theorem 4.5 gives another formula for formal 
power series corresponding to (^-distributions, namely 

(4-6) /W*tt(*0 = PmiVini 2 ))' 

where p m is given by (3.3), which turns out slightly more useful in computations in- 
volving ^iZ/i 2 than (4.4) (see Section 5). 

Corollary 4.2. If \i\ and /z 2 are tp -distributions of certain random variables, then the 
ip-moment of /iiZ/i 2 of order n e N depends on ip-moments of \i\ of orders 1 ^ k ^ n 
and p-moments of /i 2 of orders 1 ^ k ^ n — 1. 

Proof. This fact is a consequence of Proposition 4.1 and Theorem 4.5. In fact, the 
combinatorial formula for d n in the proof of Theorem 4.5 gives the explicit formula for 
N^ u (n) in terms of N^{m) and N v {k) for 1 ^ m ^ n and 1 ^ k ^ n — 1. Thus, in view 
of (1.2), a similar property holds for the (^-moments of [iLv. ■ 

Example 4.1. If e E, the lowest order first return moments A^ lZM2 (n) are given 

by 

AW(1) = ^i(l), 

N, u , 2 (2) = 7V W (2)A^ 2 (1), 

N^ 2 (3) = 7^(3)^(1) + ^(2)^(2), 

= N^)Nl{l) + 2^(3)^(2)^(1) + ^(2)^(3). 

When we return to the (^-moments /iiZ^i 2 (X n ), using (1.2), we can express them in 
terms of 'universal polynomials' Q n (fii(X), . . . , fxi(X n ), /z 2 (X), . . . , /i 2 (X n_1 )), but one 
should note that they are not homogenous as in the case of free convolutions [25]. 

5. Convolutions of measures on R + and T. 

In this section we show that the orthogonal multiplicative convolution can be defined 
for arbitrary probability measures on IR + which are not concentrated at zero, as well as 
for arbitrary probability measures on the unit circle T. 

For the following result, we refer the reader to the works of Belinschi and Bercovici 
[3,4]. 

Theorem 5.1. There is a bijection between A4r + and the class of analytic self-maps rj 
ofC\R + , such that 

(1) rj(z) = rj{z) for all z e C\R+. 

(2) \im z ^o-,z<oV(z) = 

(3) &rgrj(z) e [aigz, it) for all z e C + . 

Moreover, the map n corresponding to n e A4r + is given by the transform r/^. 
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We shall use this result to prove that formula (4.4), which was shown to hold for 
formal power series corresponding to distributions of random variables, can also be 
used to define /~tiZ/i 2 for // 2 e -Mr + and /_t 2 # So, if formal power series are replaced 
by analytic functions on C\M + . 



Proposition 5.1. If fii, /J.2 e A4r + and yU 2 is not concentrated at zero, then there exists 
a unique probability measure fi e -Mr +; such that 

(5.1) Mz) = Z ^f ] 

for z e C\R + . The measure \i will be defined to be the orthogonal multiplicative convo- 
lution of Hi and [i2, denoted fXiAfx 2 - 

Proof. For z e C\R+, let us define 

rj(z) = ^i%4^1. 

Since /_t 2 e -Mr + , transforms r/ m and 77 M2 are analytic self-maps of C\R + which satisfy 
the conditions of Theorem 5.1. This implies that r)(z) is an analytic self-map of C\M + 
since /i 2 ^ So. Condition (1) of Theorem 5.1 holds for 77 Ml and r\ m and therefore it holds 
for rj. Then, it is known that 

<c lim UtM. <c oo, 

2->0", z<0 Z 

and this limit is different from zero (it equals zero only if /i 2 = S ). Thus, we obtain 
lim 7}{z) = lim — lim V^i(v^( z )) 

2->0", 2<0 2->0",2<0 ^ (Zj 2->0",2<0 

= c- lim 77 M (77w(*)) 

2— >0 _ ,2<0 

= 

for some ^ c < go, where we used condition 2 of Theorem 5.1 for the function r]^ (if 
z — > 0~, tpn 2 (z) — > 0~ and thus i]^ 2 (z) — > 0~ and therefore Vmiv^i 2 ')) as z ^ 0~). 
This proves that condition (2) of Theorem 5.1 holds for the map 77. 

Next, for any z e C + , we have ip^z) e C + , using the integral representation of 
ipn 2 (z) given by (1.1). This implies that r]^ 2 (z) e C + , since 

„ ^ = jMf) ^(^) + KM^P; 

WJ 1 + 1U*) i + ^WI 2 ' 

and therefore 

^2 ( z ) 



^ 2 (*) 



ii+^ 2 wr 

Hence, applying Theorem 5.1 to the function rj^, we get 

argtynfawW) > argr? M2 (z), 

and thus 

argr7(z) = argz + arg77 M1 (77^ (z) ) - argr] ll2 (z) ^ argz. 
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Similarly, 



axgr](z) = argr^ (z)) + argz - argr^ 2 (z) 
= argr^ (r/ M2 (z)) - (arg?? M2 (z) - argz) 

< 7T 



since we have 



arg^(^) ^ argz and arg7/ w (r? M2 (z)) < vr, 

which follows from rj^ 2 (z) e C + . This proves that 77 satisfies condition (3) of Theorem 
5.1, which completes the proof. ■ 

For computations, it is convenient to use (4.6) and represent the transforms involved 
as continued fractions. The continued fraction representation of rj^ for compactly sup- 
ported p e A4r + is obtained from that of the Cauchy transform G^, or its reciprocal 
F^, by using the formula 

(5-2) VlM (z) = l-zF lt (^ 

If p g A4r + is associated with sequences of Jacobi parameters (a, a>), where a = (a n )n^o 
and uj = (oj n ) n ^ , which we write J ((a) = (a, u>), then, using the continued fraction 
representation of G M , we obtain 



(5.3) r]^(z) = a z + 



UJqZ 2 



OL\Z 



UJiZ 2 



a 2 z 



UJ 2 Z 2 



u 3 z 2 
a 3 z 



and a related continued fraction for p M . Conversely, this formula, together with the 
first formula of (1.7), enables us to compute the Jacobi sequences corresponding to the 
orthogonal multiplicative convolution. In the case of compactly supported measures, 
the Jacobi sequences uniquely determine the corresponding measure. 

Example 5.1. Let S a be the Dirac measure with a > and let (i e A4r + be compactly 
supported, with the associated Jacobi sequences J ((a) = (a,uj). Then p^(z) = a and 
therefore, Ps a ^n(z) = ps a (vA z )) = a > which gives 5 a Ap = 5 a . In turn, p^s a (z) = Pn(az). 
Let us denote the corresponding transformation of compactly supported measures by 
S a . In terms of reciprocal Cauchy transforms, it can be defined by a convex linear 
combination 

FsA*) = ~F Dalx (z) + (1 - -)z, 

where D\p is the dilation of measure (A by A, defined by D\p(E) = p(\~ l E). In terms 
of the transforms rj, we have r) Sa ^(z) = r] Da ^(z)/a, using (5.2). We also have r/ DafM (z) = 
rj^az). Note that S a p = (T a o D a )p, where T t -i is the so-called t-transformation of 
measures [9] . If p is compactly supported, we can write r)s afi in the form of a continued 
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fraction 



VS a u( Z ) = Q Z + 



UqCIZ 2 



uj\a 2 z 2 



OL\dZ 



Ll>2C1 2 Z 2 

1 — a 2 az 



u 3 a 2 z 2 
1 — a 3 az 

In particular, /iZ<5i = /i, thus Si is the right unit w.r.t. the operation Z on A4r + , and 

Example 5.2. If ^ = (1 — p)So + p§i and H2 = (1 — q)o~o + qb~i, then the corresponding 
reciprocal Cauchy transforms are 

z-l+p _ z-l+q 

» AZ} ~ z(z-l) ' »' AZ} ~ z(z-l) ' 

which gives 

(p — p 2 )z (q — q 2 )z 2 

^ =P+ l-(l-pK ^ =g * + l-(l-gK 
and therefore, in view of (4.6), we obtain 

(p-p 2 )qz 
P ^ {z)=P+ l + {pq-l)z 

Using (5.3), we obtain J(/i 1 Z / u 2 ) = (q,uj), with a = (p, 1 — pg, 0, 0, ...) and = 
((p — p 2 )q, 0, 0, . . .). The corresponding measure has two atoms (explicit dependence on 
p and q is rather complicated and is omitted). 

A result analogous to Theorem 5.1 also holds for the set Mj of probability measures 
on the unit circle T = {z e C : \z\ = 1}. Here, the class of self-maps rj of the open unit 
disc D = {z e C : \z\ < 1} is used, which is again related to the transforms 77^. For 
details, we refer the reader to the works of Belinschi and Bercovici [3,4]. 

Theorem 5.2. There is a bijection between Aij and the class of analytic self-maps 77 of 
D such that 77(0) = and \t](z)\ ^ \z\ for all zel. Moreover, the map 77 corresponding 
to fi is given by the transform rj^. 

Proposition 5.2. If fx 1} /j, 2 e M.f, then there exists a unique probability measure fi e 
M.f such that 

(5.4) „„ w - b^sA 

for zeD, where we understand that 77 w (0)/0 = ?7L(0). This measure n will be defined 
to be the orthogonal multiplicative convolution fii/L^. 

Proof. First of all, observe that \r)(z)\ ^ \z\ since 77^ has the same property by 
Theorem 5.3. Moreover, 77 is a quotient of analytic functions on D, and if n^ 2 (z ) = 0, 
then 
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FIGURE 3. An example of the orthogonal loop product 

thus n is a well-defined analytic self-map of O. Finally, 77(0) = • 77^(0) = 0. Thus, rj 
satisfies the conditions of Theorem 5.3 and therefore 77 = 77^ for some \x e A4j, which 
completes the proof. ■ 

Example 5.3. For any \x e M.f and a e T, we get 5 a Z/i = 5 a and /iZ<5 a = S a ti, where 
S a fj, is defined by the same equation r] Safl (z) = r] Da ^(z)/a as in the case of measures on 
K+ (on Mj, the operation D a should be interpreted as a rotation). 



6. Orthogonal loop product of graphs 

In this section we establish a relation between the orthogonal multiplicative convolu- 
tion of compactly supported probability measures and a new type of product of graphs 
called the 'orthogonal loop product'. 

As in the monotone case, we begin with recalling the definition of the orthogonal 
product of graphs [2,16]. Then we will modify this product in an appropriate manner 
to find a relation with the orthogonal multiplicative convolution. 

Definition 6.1. The orthogonal product of two rooted graphs (Qi,ei) and {§2^2) is 
the rooted graph (Q 1 h Q 2 , e) obtained by attaching a copy of Q 2 by its root e 2 to each 
vertex of Q\ but the root e±, where e is taken to be equal to e±. If its set of vertices is 
identified with V\ h V 2 := (Vf x V 2 ) u {ei x e 2 }, then e is identified with e 1 x e 2 . 

It is worth noting that the orthogonal product of graphs resembles their comb prod- 
uct. The difference is that in the comb product the second graph is glued by its root 
to all vertices of the first graph, whereas in the orthogonal product the second graph 
is glued to all vertices but the root of the first graph. 

Definition 6.2. Suppose the edges of the orthogonal product (Qi h Q 2 , e) are /-colored. 
The orthogonal loop product of rooted graphs (Qi,ei) and (^2,^2) is the rooted graph 
(Qi he <? 2 ,e) obtained from (Qi h £ 2 ,e) by attaching a loop of color 1 to all vertices 
but the root of each copy of Q 2 , and a loop of color 2 to the root of Q\. 

The justification of the above definition is similar to that in Remark 3.1. It is enough 
to deduce the glueing rules from the expressions for Ri and i? 2 , given in the proof of 
Proposition 4.1, interpreting a t as the adjacency matrix of graph Q L , tel. Also, there 
is an easy analog of Remark 3.1. Finally, it follows from [2] that the additive analog 
of the Multiplication Theorem holds for orthogonal independence. An example of the 
orthogonal loop product of graphs is given in Fig. 3 (cf. Figs. 1-2). 
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Example 6.1. In the orthogonal loop product of graphs given in Fig. 3, we have 

(Pi, w u p 2 , w 2 , fa, w 3 , /3 4 , w A ) e D 8 (e), 
where the sequence of edges of color 1 is the same as in Example 3.1, namely 

(Pi, Ih, P3, Pa) = ((e, x), (x, x'), (x', x), (x, e)), 

and the (alternating) f-walks attached to the vertices that order) are: 

wi = w 3 = (x,x) 2 , w 2 = (x',x') 2 and u>4 = (e, e) Note that this d-walk has lenght 8. 
Again, as in the case of d- walks on the comb loop product, if any of the loops at x or x' 
is replaced by an alternating f-walk, we shall still get a rooted alternating d-walk. Of 
course, the simplest rooted alternating d-walk is given by (Pi, wi, P4, w±) e D\(e) and 
(/3 ,w 4 ) e D 2 (e), where (3 = (e,e) l . 

Theorem 6.1. The Multiplication Theorem holds for orthogonal independence, the as- 
sociated loop product Gi \- e G 2 , an d the multiplicative convolution [i\L[i 2 . 

Proof. The proof is similar to that of Theorem 3.1, but the combinatorics is based 
on a different formula for formal power series, namely 

which leads to the combinatorial formula 

n 

JW(n) = 2 N^r) J N^N^h) . . . N^k^), 

r=l ki+k2 + ...+k r -i=n— 1 

where it is assumed that hi, k 2 , . . . , k r _x are positive integers. Recall that in the or- 
thogonal product Qi h Q 2 there is one copy of Q\ (with e\ identified with the root e of 
the product graph) with a copy of Q 2 attached by its root to every vertex x of Q\ but 
the root e\. Therefore, each d-walk w e D 2n (e) is a sequence of r edges of Gi, which 
themselves must form an f-walk c = (e, vi, v 2 , . . . , iV-i, e) of color 1, interlaced with 
alternating f-walks Wi e F(vi), 1 ^ i ^ r — 1 and a loop of color 2, w r = (e, e) 2 . Note 
that this is the only way to produce a rooted alternating d-walk since the only edge of 
color 2 incident on e is the loop and therefore, in order to get an alternating walk, the 
first f-walk which begins and ends with an edge of color 1 must be followed by the loop 
at e of color 2 to make a double return to e. The contribution from each product of 
type 

N fll (r)N, 2 (k 1 )N, 2 (k 2 )...N, 2 (k r _ 1 ) 

to the RHS of the above formula is equal to the number of all such f-walks c of color 
1 which have r edges and are interlaced with r — 1 f-walks of lenghts 2k\ — l,2k 2 — 
1, . . . , 2A; r _i — 1 attached to all vertices of c but one (we choose this vertex to be the root 
since we want the considered walk to be an f-walk). The summation over all 1 ^ r ^ n 
and 2ki + 2k 2 + . . . + 2A; r _i — r + 1 = 2n — 1 indicates that the first f-walk ui in the 
d-walk w = (ui, u 2 ) is of lenght 2n — 1, and u 2 is the loop of lenght 1. The summation 
over 1 ^ r ^ n gives exactly the cardinality of D 2n (e), which finishes the proof. ■ 
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Example 6.2. Let us apply the formulas of Example 4.1 to the enumeration of rooted 
alternating d-walks on the graph Q\ he G2 in Fig.3 (we keep the notation of Example 3.2 
for spectral distributions /ii and \x-i and thus we get the same values of the (fc)'s and 
the N^^jys). Using Example 4.1 and Theorem 6.1, we get ^(e) = D±(e) = D§(e) = 
Ds(e) = 1, which can be verified directly by counting rooted alternating d-walks on 
Q\ he &2- 



7. Subordination operators 

In analogy to the additive case [16], where we introduced and studied operators 
related to the subordination property for the free additive convolution, we shall now 
present an analogous approach to the subordination property for the free multiplicative 
convolution. We will mainly refer to sets and Aij, which correspond to positive 

and unitary operators, but the operatorial subordination will also hold for all bounded 
operators. 

Let (W t ,£ 4 ), where 1 e /, be Hilbert spaces with distinguished unit vectors. Then 
their Hilbert space free product (7ii,£i) * (7^2, £2) is (W, where 

00 

(7.1) h = C £0© © K®Hl®...®Hl, 

with = Tt L © C£ t and £ denoting a unit vector (canonical scalar product is used). 
For any h e H c , denote by h° the orthogonal projection of h onto H° L . Moreover, let 

(7.2) WW( t )= © K®<®---®K> 

(7.3) /CW(0= © K®Hl®...®Hl, 

for l e I and n e N, and, for convenience, we set H^(l) = IC^(l) = C£ with the 
canonical projection P$ : — > C£. 

Since our index set I consists of two elements, the above notation gives identifications 

nodd, 
JC^^T) n even. 



(7.4) H {n \i) 



Nevertheless, it is convenient to use both sequences, (7.2) and (7.3), as well as direct 
sums 

00 00 

(7.5) W(t) = ®W (n) (0 and = ©/C^(i), 

where tel. 

As in our previous work, we decompose the free product of Hilbert spaces as (two 
different) orthogonal direct sums 

00 

(7.6) H = ®n {n - l \i)®l-L {n ^), 

71=1 
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where i e I, and denote by 

(7.7) P t (n) : H — (i) ® H (n) (l) 

the associated canonical projections. Finally, we define the so-called vacuum state 
¥>(■)=<■& °n B(W). 

Of particular interest will be s-/ree products of Hilbert spaces, (7i t ,£ t ) and (7i r , £ t ), 
denoted (H L ®Hj, £) and defined [16] as the pair (X t ,£) where 

(7.8) /c t = cee/c(i), 

and tel. The direct sum decompositions 

(7.9) /c t = © w^^weww^ 

n odd 

hold for each tel. 

Let x e B(Hi) and y e Bi^K-i) be fixed random variables. The corresponding free 
random variables X\ := \(x) and X 2 := A(y) are elements of B(7i), where A is the free 
product representation on H, and can be decomposed according to (7.6) as 

oo 

(7.10) X t = £ 

3 = 1 

where X L (j) = P L (j)X L P L (j), j e N, can be viewed as replicas of x and y, respectively, 
where l e L Using (7.6) and (7.9), we also have the decompositions 

oo 

(7.H) 1 = 2^0') and l Ki = P.U) 

3=1 j odd 

of the units in B{H) and S(/C t ), where te/, respectively. Here, and also in the sequel, 
we denote by the same symbols, X t (j) and P L (j), the corresponding operators on /C t . 

We shall use representations of free random variables as 'orthogonal series' with the 
unit singled out, namely 

oo 

(7.12) X L = 1 + J> t (j) 

3=1 

where x L (j) = X L (j) — P t (j) and i e /, j e N. If no confusion arises, we will distinguish 
the first term in the above series by a special notation and write x t = x t (l). Let us 
remark that the above form is suitable for the study of multiplicative convolutions, 
where units play a special role (roughly speaking, they have to be subtracted when it 
comes to proving some kind of independence). 

We begin with the decomposition of the product X-X b which corresponds to Eq.(1.3) 
for / ui,//2 e -M-R+, and therefore involves positive operators. 

Theorem 7.1. Let X L , i e I, be positive random variables which are free with respect 
to ip, and denote t L = 1 + x L , T T = \IX L — x l X T \/X i — x L . Then X T X L has the same 
ip- distribution as ^Jt i T T ^Jt L and the pair (t L — 1,T T — 1) is monotone independent with 
respect to (p. 
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Proof. Note that T T introduced above is well-defined since 

oo 

X t -x t = P t (l) + 2^0') 

i=2 

is positive as an orthogonal direct sum of positive random variables. Using a similar 
argument, we get positivity of 

oo 

t t = x t (i) + 2*Uj)- 

J =2 

Computations of square roots give 



V^(1) + 2^(J) and v / ^ = P t (l) + ^ v ^(]j 1 

where sJX L (j) = P L (j)\{sfz)P L {j) for j e N, with z = i or z = ?/ if t = 1 or t = 2, 
respectively, which leads to 



Now, we use the fact that X-X t has the same (^-distribution as ^/X~ L X-^/X~ i , which 
completes the proof of the first part of the theorem. 

We will show now that the pair (t L — 1,T T — 1) is monotone independent w.r.t. ip. For 
simplicity, we denote x = t b — 1 and y = T- — 1. 

Case 1. If y n is in the middle of the moment, we compute ip(wixy n xw 2 ), where neN 
and wi, w 2 e alg(x, y). Note that the range of x is C 7i° L , therefore, we only need to 
find the action of y n onto £ and h e HP L . Using the explicit form of \JX L — x~ t , we get 

oo oo oo 

V n i = ((1 + £ r*)(l + £*j0(1 + 2 r*)-l) B * 

k=2 j=l k=2 

= V (y n )t mod/C(0, 

where r k = ^JX L {k) — P L {k) and z 3 - = X T (j) — P T (j) and therefore r fc £ = for any k ^ 2. 
Similarly, we get 

y n h = (y n h, h)h mod /C(l) /C«(i) 

for any h e of norm || /i ||= 1. Note that /i plays the role of a cyclic vector for the 
restriction of the free product representation 7Ti * n 2 to the unital algebra generated by 
y. Thus (y n h, h) = tp(y n )h, which gives 

(p(w 1 xy n xw 2 ) = ip(y n )ip(w 1 x 2 w 2 ), 

i.e. the required condition for monotone independence. 
Case 2. If y n is at the end of the moment, we get 

(fiwixy") = (p(y n )(p(w!x), 

using the relation for y n £ obtained above and the fact that K,(i) a Kerx, which finishes 
the proof. ■ 

Corollary 7.1. //, in addition, X L has tp- distribution ji L # 5q for l e I, then the 
p>- distribution ofT L is given by the s-free multiplicative convolution fi L \Alh- 



OPERATORS RELATED TO SUBORDINATION 



29 



Proof. First, note that the variable t b has (/^-distribution \i L . In view of (1.5) and 
Theorem 7.1, we have 77 w g M2 = ° ^ = ° ^oi, where cji and o" 2 are the </?- 
distributions of Ti and T 2 , respectively. From the definition of the s-free multiplicative 
convolution, it follows that o~i = 0yU 2 and a 2 = /i 2 \Z\fJ>i, which proves our assertion. 
■ 

Remark 7.1. In analogy to the terminology used in analytic subordination [8,23], 
we can say that X-X L is 'subordinate to t L and with T- and T L , respectively, being 
the corresponding 'subordination operators'. Let us remark, however, that operators 
T L) i e J, are not good candidates for 'subordination branches' of the product of free 
random variables since they are not suitable for further decomposition. In that sense, 
they are not multiplicative analogues of the 'additive subordination branches' [16], or 
the branches of the free product of graphs introduced by Quenell [20]. In the next section 
we will introduce and study such analogues, whose decompositions will correspond to 
decompositions of s-free multiplicative convolutions. 

A similar operatorial subordination result can be established for bounded operators 
which includes unitary operators (the latter are related to A4j). 

Theorem 7.2. Let X L , eel, be bounded random variables which are free with respect 
to (p. Then, for each eel, the operator X-X i has the same <p -distribution as Q T q L , 
where q L = 1 + x v and Q- = X-(X L — x b ) . Moreover, the pair (q L — 1, Q- — 1) is monotone 
independent with respect to ip. Finally, if X L , t e I, are unitary, then the operators 
q t , Q b , i e I, are unitary. 

Proof. The proof of the statements concerning bounded operators is similar to that 
of Theorem 7.1. Therefore, we shall just prove unitarity of Q L and q t . Writing operators 
q L and X L — x L in the form of orthogonal series 

oo oo 

q L = X t (l) + J] P L (j), X t -x t = P t (l) + J] X t (j), 

we obtain 

oo 

q iq : = x L (i)x;(i) + Y j P.(j) = i 

3=2 

and 

oo 

(x L - x L ){x b - Xb y = p L {i) + J x L (j)x:(j) = i. 

3=2 

Similarly, q*q L = 1 and (X L — x i )*{X i — x b ) = 1. This proves unitarity of q L and X L — x L , 
from which we obtain unitarity of Q L . The remaining part of the proof is similar to 
that of Theorem 7.1. ■ 

Corollary 7.2. Let /i t and a L be the <p -distributions of X L and Q L , respectively, for 
tel. Then the ip -distribution of X-X b is given by fi L Q a- for each tel. In particular, 
if the operators X L are unitary, then o L = /i t for each tel. 

Proof. The first statement follows from Theorem 7.3 and the definition of the mono- 
tone multiplicative convolution of distributions. The proof of the second statement is 
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similar to that of Corollary 7.2 (the s-free convolution of measures from M* is used). 



8. Subordination branches 

In order to introduce operator- valued 'multiplicative subordination branches' asso- 
ciated with the product of bounded free random variables, let us introduce random 
variables 

(8.1) R t (m) = l+ x t (j) 

jeJ(m) 

where the notation J(m) = {m, m + 2, m + 4, . . .} is used for index sets and, by abuse of 
notation, 1 denotes 1^ if m is odd, or if m is even. Note that operators R L (m) are 
elements of B{KL L ) or B(JC T ), depending on whether we have m odd or even, respectively. 

Therefore, on each of the 's-free Fock spaces', /C t , i e I, we get an 'interlaced' sequence 
of operators 

(8.2) R t (l), R T (2), R t (3), R-(A), ...e B{K, b ). 

Each of these two sequences is used in the definition of one sequence of operator-valued 
'multiplicative subordination branches'. 

Proposition 8.1. If X±,X 2 are positive (unitary), then operators R L {m), where Lei 
and meN, are positive (unitary). 

Proof. This fact can be easily checked by decomposing the units according to (7.11). 
We choose to show unitarity of 

r = X^) and r = ^+ S^) 

j odd j even 

for unitary X t . Since X t {j)X*(j) = P L (j)X L X*P L (j) = P L (j) for each te/andje N, 
we get 

R: dd (R: d r = 2 x^wv) = 2 = ^ 

j odd j odd 

Rr{Rrr = p € + 2 x L (j)x;(j) = p € + 2 p t (j) = i. 

j even j even 

Unitarity of R t (m) for arbitrary m can be proved in similar way. Positivity of R L {m) 
follows easily from the appropriate decomposition of the unit. ■ 

Definition 8.1. Let X t , i e I, be random variables which are free w.r.t. if. Random 
variables 

(8.3) U,(m) = R t (m)R T (m + l), 

where t e I and m e N, will be called multiplicative bounded subordination branches. 
If X L , Lei are unitary, then they will be called multiplicative unitary subordination 
branches. 
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In addition to branches U b {m), we shall need operators obtained from U L [m) by 'order 
reversal'. Namely, let 

(8.4) V t {m)=R T {m + l)R b {m) 

where Lei and m eN. Again, we use a simpler notation for m = 1, namely V L = V t {m). 

Definition 8.2. Let X L) i e I be positive random variables which are free w.r.t. ip. 
Random variables 

(8.5) Y L {m) = y/R T {m + l)R L {m)y/ R T {m + 1), 

where t e I and meN, will be called multiplicative positive subordination branches. 

For simplicity, we will also use the terms: 'bounded branches', 'unitary branches' and 
'positive branches'. As in the case of operators R L {m), where m e N and t e 7, we get 
two sequences of alternating (bounded, unitary) branches 

(8.6) U L (1), U T (2), U L (3), U T (4), ... e B(K, U ) 

for each Lei, and similar sequences of positive branches. Recall that existence of 
'interlaced branches' was also observed for the additive (self-adjoint) subordination 
branches [16]. 

Of particular interest are branches of 1st order, and it is therefore of advantage to 
use a special notation for the first two operators in each 'interlaced' sequence given by 
(8.1). For l e I, we shall use 

(8.7) Rf d = 1+Y i x l (j) and R- = 1 + £ x L (j) 

j odd j even 

where again, the unit 1 has to be interpreted as ljc t and l/e r , respectively. Note that 
although we have two different units here, the operators which appear in the same 
'interlaced' sequence contain the same unit. 

Remark 8.1. Let us recall that in the additive case [15, Theorem 8.4], the branches 
B\ and B 2 of the sum X± + X 2 can be decomposed as 

(8.8) B 1 = S{ dd + S? cn , 

where Sf d = £ jodd Xi(j), S? cn = 2 jeven X 2 (j), and that the pair (Sf^S^) is s-free 
w.r.t. the pair of states (ip, ip), where ip is associated with any unit vector £ e H®. Note 
that we have Rf d = S'J ,dd and i?2 vcn = P% + S^"" 1 , but we prefer to have a new notation 
to maintain a uniform style for 'subtracting units'. 

This, in turn, leads to a special notation for branches of 1st order, U t = U t (l), 
Y L = Y L (1), Lei. For instance, 

Y l = Y l (l) = y/R^Rl^y/R^, 

Y 2 = y 2 (l) = ^/RY cn R° 2 dd ^/Rr n - 

Many computations can be reduced to branches of 1st order since their (^-distributions 
agree with the distributions of branches of higher orders with respect to suitably chosen 
states. 

Let us examine the (/^-distributions of branches of first order. In fact, we will show 
that Y L and U L have the same ^-distributions as T L for given tel. 
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Lemma 8.1. The variables Q L , U L and V L have the same if -distributions for any given 
tel. If X L , i e I, are positive, then the variables T L ,Y L ,U L and V L have the same 
ip- distributions for any given tel. 

Proof. We shall give the proof in the more difficult case of positive operators (the 
proof for the case of bounded operators is similar). For n ^ 1, we have 

<P(YT) 



since ^JR^ CU ^ = £ (the only summand in the direct sum decomposition of this square 
root which gives a non-zero contribution is 1-(1))- This shows that Y t and U t have the 
same (/^-distributions. 

Let us now compare the moments of T L with those of Y L . We have 

V9(77) = <p (V*r - x T (X L (X T - x T )) n - l X^X T - x r ) 
= v ((X L (X T -x T )) n ) 

since ^JX- — x-C, = £. To fix attention, suppose that i = 2 and denote, for convenience, 
Zk = xi{k) and y k = X2(k). Now, let us observe that in the finite set of variables of type 
Zk and yj, where k ^ 2 and j ^ 1, which are used when calculating X 2 {X\ — xi) n-1 X 2 £ 
there are no z^'s for odd fc's, or y/s for even j's. This is because the only variable 
which may give a non-zero contribution when acting on £ is yi, then only y x and z 2 
may follow, of which the first operator can be followed by y\ or z 2 , whereas the second 
- by z 2 , Hi or Vz-> etc. This leads to 

<p(p(x 2 , x 1 - Xl )) = <p(P(R? d , R e D) 

for any polynomial P in two noncommuting variables. This shows that T 2 and Y 2 have 
the same (^-distributions. 

Finally, we will show that U 2 and V 2 have the same (^-distributions. Using the 
notation used in the previous paragraph, we have 

tpm-ir) = ^(2 z k + Y,Vj+ 2 

k even j odd k even, j odd 

<p ((v 2 - 1)") = </>« I> + 2 + 2 

fc even j odd fc even, j odd 

b-fc|=i 

for any natural n, and therefore, we need to show that the right-hand-sides of the 
above equations are equal to each other. Note that the products of variables z k and yj 
in words w which appear in mixed moments <p{w) giving non-zero contributions to the 
above moments satisfy the following three conditions: 

(i) whenever z k stands next to z m or y k next to y m , it holds that k = m, 

(ii) whenever yj stands next to z k , it holds that \j — k\ = 1, 



= (y/Rr* (R° ddR 7 cn ) n ~ lR T d \/PF cn ) 

= (p ((Rl dd Rj cn ) n - l R° dd ) 
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(iii) the word w begins and ends with the letter y 1 . 
Let us denote by W the set of words in letters 

L = {z k , Dj : k even, j odd} 

subject to conditions (i)-(iii). In turn, let W and W", respectively, denote the sets of 
all words in letters 

V = {z k ,y j ,z k y k - 1 ,z k y k+1 : k even, j odd}, 
L" = {z k ,y^y k -xz k ,y k+x z k : k even, j odd} 

subject to conditions (i) and (iii), where products of type z k yj and yjz k are treated 
as letters. When they are used to form words, they are denoted by [z k yj] and [t/j-Zfc], 
respectively. Note that W a W and W a W" . More importantly, there is a bijection 

r : W -> W" 

defined by 'shifting brackets to the left'. Namely, r is uniquely defined by r(u>i) = W\ 
and the recursion 

T{w x y m z k \z k yf\w 2 ) = w 1 [y rn z k \z p k T(y j w 2 ) 

for any w\ e W, w 2 e W, p ^ and \j — k\ = 1, \k — m\ = 1, i.e. we assume that [z k yj] 
is the first such pair (counting from the left). In this recursion, we shift the brackets to 
the left from this pair to the closest possible pair y m , z k lying to its left (in view of (iii), 
such a pair must exist). It is not hard to see that its inverse is given by t _1 (m 2 ) = u 2 
and the recursion 

^ X {w'x\V3Z k \zly m w' 2 ) = ^{w^y^zllzkyr^w'z 

where w[ e W", w' 2 e W, p ^ and \ j — k\ = 1, |fc — m\ = 1. Again, note that for every 
such pair [yjZ k ] there must exist y m standing to its right in view of (iii). Clearly, for 
every w e W it holds that (p(r(w)) = <f(w). This completes the proof. ■ 

We are ready to prove a theorem, which will lead to a relation between the distribu- 
tions of two consecutive branches. For that purpose, introduce operators 

(8.9) Z t (m) = l + x L {m) 

with the same convention concerning the units as before (1 = 1^ if m is odd and 
1 = l)c- if m is even), as well as operators obtained from the subordination branches 
by subtracting the (appropriate) unit: 

(8.10) y L (m) = Y L (m) - 1 
for any meN and tel. 

Theorem 8.1. For given m e N and i e I, let (p,ip be states associated with any unit 
vectors ( e H^ 171 ^^) and(' e H^ m \l). Then, operators U L {m) and Z t (m)U T (m + l) have 
the same (p -distributions. Similarly, if X b , i e I, are positive, then operators Y L (m) and 
Z L (m)Y T (m + 1) have the same ip -distributions. Moreover, the pairs {x L {rn),Uj{m + 1)) 
and (x t (m), yxijn + 1)) are orthogonal with respect to (tp, if)) for each i e I and m eN. 
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Proof. We shall prove the more difficult case of positive random variables. Without 
loss of generality, let i = 1 and, for notational simplicity, consider the case m = 1 (the 
proof of the general case is similar). Denote 

Z 1 = Z 1 (l), R l = R^l), R 3 = ^(3), R 2 = R 2 (2), Y 2 = Y 2 (2). 

In the first part of the theorem we need to prove that Y\ has the same ^-distribution 
as Z\Y 2 . Using orthogonal decompositions and 'square root calculus', we obtain 

Z 1 =X 1 (1)+ J] Pi (j), and y/R 3 = P 1 (l)+ £ VMt) 



ieJ(3) 



ieJ(3) 



which gives 

and therefore 
Therefore, 



j odd 



Y 1 



P2P1V R-2 



R 2 y RzZxJ R3V R-2- 



V {WR2 



R 3 Z 1 ^/R 3 ^/R 2 
R , i{Z\Y 2 ) n X Z\\J~R 3 



n-l 



Ziy R3' 
R2) 



'Ri 



since Y 2 ^ = £ and also 



v(\fR~2' 
ip ((z 1 Y 2 y- 1 z 1 ) 

V ((ZiY 2 ) n ) 

^RliVRii = C Therefore, the (^-distribution of Y\ agrees with 
the (^-distribution of Z{Y 2 . 

In the proof of orthogonality, denote, as before, x\ = £i(l), and y 2 = Y 2 — 1^. Using 
orthogonal decompositions and 'square root calculus', we can write 

Y 2 = y / P^P 2 y / i?3 

= (A(i)+ 2 VmT))(p,+ X x 2 (k))(p 1 (i) + X VMf)) 

jeJ(3) keJ(2) jeJ(3) 

and therefore Y 2 £ = £, which implies that y 2 ^ = 0. Therefore, we get (p(wy 2 ) = and 
therefore, by taking the adjoints, ip(y 2 w) = for any w e alg(xi,y 2 ). This gives the 
first orthogonality condition. Now, let us prove the second orthogonality condition, i.e. 

ip(w l x l ylxiw 2 ) = i)(y 2 )((p(wix\w 2 ) - w{w 1 x 1 )y{x 1 w 2 )) 

for any k e N and Wi,w 2 e alg(xi, y 2 ). Note that the 'lowest order term' in the 
expression for y 2 , which turns out to be of special importance when x\ stands next to 
y 2 , is of the form 

P 1 (1)X 2 (2)P 1 (1) = P 1 (1)P 2 (2)X 2 P 2 (2)P 1 (1). 

Recall that P 2 (2) is the projection onto Ti\ © (TL 2 ®T~Li) and Pi(l) is the projection 
onto C£©ft?, thus Px(l).P 2 (2) = P n o, and therefore we get 

XlVT x l = XlPnoy^P-HOXi, 



OPERATORS RELATED TO SUBORDINATION 35 

for m 1. This gives 

ip(wixiy™xiw 2 ) = y(w 1 xiP H ox 1 w 2 )il){y™), 
since, for any £ e 7i\ of norm one, we have 

y?( = V>(y 2 m )C mod/c(2)e/c«(2), 

where if) is the state associated with Note that -0 gives the same moments of the 
variable y 2 irrespective which C is taken and they agree with the corresponding moments 
in the state ip 2 . Finally, since P n o = P c ^ H o — P c ^, we have 

ip(w 1 x l P n ox 1 w 2 ) = ip(w 1 xlw 2 ) - ip(w 1 xi)ip(x 1 w 2 ), 
which completes the proof of the second orthogonality condition. ■ 

Corollary 8.1. Let /x t and a L be ip -distributions of (bounded operators) X L and U L , 
respectively, for each l e I, and let V>m-i be the state associated with any unit vector 
( e W (m_1) (t), where m e N. Then 

(1) the ip m -i -distribution of U L (m) agrees with a L for every m ^ 2 and eel, 

(2) it holds that o L = fi t Aa T for each i e I, 

(3) if Xi,X 2 are positive and /i 2 # S , then a L agrees with the s-free convolution 

\Z\^t of probability measures on R + for each i e I, 

(4) if X 1 ,X 2 are unitary, then a L agrees with the s-free convolution ii l \Z\ of prob- 
ability measures on T for each tel. 

Proof. Assertion (1) follows from the definition of operators U L (m). Then, (2) is a 
consequence of (1) and Theorem 8.3. Finally, it follows from Theorems 7.3 and 8.3 that 
the (^-distributions of U L satisfy the subordination equations (2.1). Uniqueness of the 
subordination functions for the cases considered in (3) and (4) give the assertions. ■ 

Corollary 8.1, together with Corollary 4.6, naturally lead to the definition of a se- 
quence of iterations of the s-free multiplicative convolution. Thus, for given /ii,/i 2 6 S 
let 

(8.11) /UiZi£t 2 = ^1^2, /il^nM2 = /ilZ(/i 2 Z n _i//i), 

for n ^ 2 (of course, this sequence can also be defined by for any measures, for which 
the operation Z has been defined). 

Corollary 8.2. If fii, fi 2 e A4r + \{Sq} are compactly supported, then 

(8.12) w - lim (AtiZ n /x 2 ) = A*i A*2 

n— >oo 

(8.13) w - lim (fi! Q (/i 2 Z n /zi)) = a*i fJ, 2 , 

n— >oo 

and the corresponding r\ -transforms converge uniformly on the compact subsets of C\R + 
to the n-transforms of the limit measures. An analogous result holds for /i 2 e 

Proof. In view of Corollary 4.6, for fixed m e N, (yUiZ„yU 2 )(m) = (/iiZ m /i 2 )(m) for 
all n ^ m, from which we obtain convergence of moments as — > oo for any distributions 
fj,i, /i 2 g S. Since the measures are compactly supported, this implies weak convergence 
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of the corresponding measures. Therefore, the 77-transforms converge uniformly on 
compact subsets of C\R+. Moreover, by Corollary 8.1, the limit distributions agree 
with the distributions of the positive subordination branch Y]_, which is p\ 0// 2 . The 
arguments for the weak convergence of the second sequence are similar. An analogous 
proof holds for p\ , pi £ * . I 

Remark 8.2. Informally, the weak limits of Corollary 8.2 can be written in the following 
form: 

Pi A*2 = /iiZ(/i 2 Z(/iiZ(/i 2 Z(. . .)))), 

Hi fl2 = /!! O (/i 2 Z(//iZ(/i 2 Z(. . .)))), 

whereas their transforms in the 'continued composition form': 

P/n0^O) = Pn( z Pv*{ z Pn( z Pvu{- ..)))), 
{zp^{zp^{zp m {. ..)))), 

where the right-hand sides are understood as the uniform limits on compact subsets of 
C\R+ or D. Actually, these formulas can be used to compute some examples, including 
simple examples of p\ \x\ p 2 (without using S-transforms) . 

Remark 8.3. By Corollary 8.2, one can define the s-free and free convolutions of any 
compactly supported measures p±, p 2 on C as the weak limits of the form (8. 12)- (8. 13). 
This allows us to use operatorial subordination results of Sections 7-8, including those 
which concern distributions of products of free bounded random variables as well as 
distributions of bounded branches and express them in terms of s-free convolutions, 
denoted with the same symbol and understood as weak limits of type (8.12)-(8.13). 
Using the weak limits, we also get 

(8.14) 5 \Z\p = 5o and p 8 = S^ X ) 
and 

(8.15) p 5 = 5 = 5 /i, 

for any compactly supported probability measure on C, where we used the well-known 
relations p, Q 5 = 5 = 5 Q S^ X )- 

Example 8.1. Using Corollary 8.2 and the results of Example 5.1, we obtain p\Z\5 a = 
pA5 a = S a (p) and 5 a p = 5 a Ap = 5 a for compactly supported p e A4r + and a > 0. 
This gives 

S a \x\p = 5 a O S a p = D a p, 
in view of (1.5), since the corresponding ^-transform is of the form 

Vs a (vs a „(z)) = VD a p(z), 
where we used r)s a (z) = az and rfs afl (z) = VD a n(z)/a. On the other hand, 

P $a = A* O $a = D a p, 

since f]^ a (z) = rj^az) = r] DafM (z). In particular, we have S a 5b = S a for a, b > 0. 
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Example 8.2. Using the results of Example 5.3, we obtain 



A* 5a = S a jjL and 5 a A* = ^< 



for compactly supported /x e A^t and a e T, which gives 



jtt@5a = D a /i = 5 a \x\/i, 



by the same arguments as those used in Example 8.1. 

Let us finally show that the multiplicative subordination branches are related to the 
notion of s-free independence, as it was the case of their additive counterparts. For that 
purpose, let us recall this concept [16]. 

Definition 8.3. Let (A, <p, ip) be a unital algebra with a pair of linear normalized 
functionals. Let A\ be a unital subalgebra of A and let A 2 be a non-unital subalgebra 
with an 'internal' unit 1 2 , i.e. l 2 b = b = bl 2 for every b e A 2 . We say that the pair 
(Ai,A 2 ) is free with subordination, or simply s-free, with respect to ((p,ip) if ip(l 2 ) = 1 
and it holds that 

(i) Lp(ciia 2 . . . a n ) = whenever a,j e A® and %\ # i 2 ^= ■ ■ ■ ^ i n 

(ii) (p(wil 2 w 2 ) = (f(wiw 2 ) - ip(w 1 )ip(w 2 ) for any w 1 ,w 2 e dlg(Ai,A 2 ), 

where A® = A\ n kery? and A 2 = A 2 n ker^. We say that the pair (a, b) of random 
variables from A is s-free with respect to (ip, if)) if there exists 1 2 e A such that (Ai,A 2 ), 
where Ai := alg(l, a) and A 2 := alg(l 2 , b), is s-free w.r.t. (</?, -0). 

Proposition 8.2. T/ie pazr (i?° dd — 1, B^ cn — 1) is s-free with respect to the pair (<p, ip), 
where if) is the state associated with any unit vector ( e HP L , and the projection associated 
with the second variable is 1 2 = 1 — P^. 

Proof. For simplicity, choose t = 1 and denote r 1 = Rl dd — 1, r 2 = R C 2 Cn — 1. We have 



Clearly, the definition of 1 2 immediately gives ip(l 2 ) = 1 as well as condition (ii) of 
s-freeness. Let us show that condition (i) also holds. Let A± = alg(r l5 1) and A 2 = 
alg(r 2 , 1 2 ). If in (i) we take a n e A 2 nker-0, then A(a n )£ = 0, where A is the free product 
representation, and hence (i) holds. If a ra ,a n _ 2 , . . . e Ai n kenp and a n _i,a„_3, . . . e 
A 2 n ker-0, then the GNS representation space for such a mixed moment is /Q (none 
of a n _i, a n _3, . . . ever gets to act on £). Therefore, in the computation of this kind 
of moment, one can replace r\,r 2 and 1 2 by X X ,X 2 and 1, respectively. But then (i) 
follows from the usual freeness condition. ■ 




and 




j odd 



j even 



9. Free product and s-free loop product of graphs 



In this section we prove the Multiplication Theorem for s-free independence and free 
independence. 

Let us recall the definition of the s-free product of rooted graphs [2,16]. 
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Figure 4. An example of the s-free loop product 

Definition 9.1. By the s-free product of rooted graphs (Qi,ei) and (G2,e 2 ), denoted 
(G1OG1, ei), we understand the inductive limit of the sequence (£>i(m), ei) m( =N of rooted 
graphs, where B±(l) = Gi |— G2 and B b (m) is obtained from B L {m — 1) by attaching by 
its root a copy of Gi (if m is even), or a copy of G2 (if m is odd) to every vertex of the 
difference B L (m — l)\B L (m — 2). 

Remark 9.1. It can be seen that (Bi(m), ei) m( =N is a sequence of growing graphs with 
the root kept to be e L , and therefore, the inductive limit exists. In a similar way we 
define (G 2 © Gi, e 2 ). It is easy to see that 

(9.1) B^m) =Gih B 2 (m - 1) and B 2 (m) = G2 h B x {m - 1), 

which can also serve an inductive definition of the considered sequences. These for- 
mulas are important since they correspond to the decomposition of the s-free additive 
convolution of the form 

(9.2) m ffl /i 2 = m h (fi 2 \- (hi h (fi 2 h . . .))), 

where fi L is the spectral distribution of Q L [2]. 

Using the decomposition (9.2), we proved the Addition Theorem for s-free indepen- 
dence in [2]. Our first goal in this Section is to prove the corresponding Multiplication 
Theorem. For that purpose, we introduce a new notion of graph product, called the 
's-free loop product'. We follow it up with a proposition which justifies the definition. 

Definition 9.2. Suppose that the s-free product of graphs (Gi ©<?2,ei) is naturally 
colored. The s-free loop product of (Gi, e\) and (G2, e 2 ) is the graph (£i(B^2, ei) obtained 
from (Gi © G2, ei) by attaching a loop of color 2 to the root e\. 

Example 9.1. Consider the s-free loop product of graphs given in Fig.4, which can be 
viewed as a 'binary tree with loops' (the product graph is of course infinite, only our 
picture is truncated at distance 5 from the root). Note that graph Gi has a loop at the 
root, and so do copies of Gi in the product graph, whereas G2 does not have any loops. 
We produce Gi (BG2 in the inductive way (Definition 9.1 is used). Now, in order to 
obtain the s-free loop product, we attach one additional loop of color 2 to the root (we 
draw this loop larger than loops of color 1). On the level of the adjacency matrices, this 
additional loop corresponds to the projection as the proof of Proposition 9.1 shows. 
Note that the main purpose of attaching this extra loop is to fit the s-free multiplicative 
convolution into the general scheme of the Multiplication Theorem, where counting 
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rooted alternating d-walks leads to the moments of multiplicative convolutions. Direct 
computations give, for instance, D 2 (e±) = 1, D 4 (ei) = 0, D 6 (ei) = 4, D 8 (ei) = 0. 

Proposition 9.1. Let A L be the adjacency matrix of Q L , where Lei. The adjacency 
matrix of the s-free loop product of Q 1 and Q 2 takes the form 

(9.3) A(G 1 (Be Q 2 ) =Ri + R 2 

where 

(9.4) #i = 2 A ^ and R * = P « + X MQ- 

j odd k even 

Moreover, the ip e - distribution of Ri is the if) -distribution of R 2 is fi 2 and the pair 
(Ri — 1, R 2 — 1) is s-free independent w.r.t. (<p e , ip), where if) is the state associated with 
any unit vector ( e 7i\ = l 2 (V®). 

Proof. It follows from [2,16] that the adjacency matrix of the s-free product of graphs 
is of the form 

A(g 1 ®g 2 ) = s 1 + s 2 , 

where 

S, = A 1 (j) and S 2 = J] A 2 (k). 

j odd k even 

By Definition 9.1, the adjacency matrix of the corresponding s-free loop product is 
obtained by adding the projection P to the adjacency matrix of the subgraph of color 
2, namely S 2 , which corresponds to attaching the loop of color 2 to the root. This gives 

A(g 1 ® e G2) = R1 + R2, 

where 

#1 = Si and R 2 = P 6 + S 2 . 

In view of Proposition 8.2, the pair (i?x — 1, R 2 — 1) is s-free independent w.r.t. the pair 
((p,ip), where ip is the state associated with any unit vector ( g (here, 1 denotes 
the unit on l 2 (V), where V is the set of vertices of £1 g 2 , which corresponds to K-i). 
Moreover, it is easy to check that the i/vdistribution of R\ is fi\ and the ^-distribution 
of R 2 is n 2 . This completes the proof. ■ 

Thus, the first part of the Multiplication Theorem for s-free independence is proved. 
We now need to find a connection between the 'first return moments' of R 2 R\ and 
cardinalities of the sets D 2n . In that context, note that in order that the set D of all 
rooted alternating d-walks on g x ©£ g 2 be non-empty, either £1 or g 2 must have a loop 
at ei or e 2 , respectively. In fact, otherwise, starting an alternating walk from the root 
of g± ®g 2 , one cannot return to it since the m-th edge of that walk must be of different 
color than the m — 1-th edge and thus it must belong to Bj(m)\Bj(m — l), and therefore 
its distance from the root (understood as the distance of the vertex closer to the root) 
equals m — 1 and thus tends to 00 as m — > 00. 

Theorem 9.1. The Multiplication Theorem holds for s-free independence, the associ- 
ated loop product £1 (Bi g 2 , and the multiplicative convolution fx± 0^2- 
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Proof. For notational simplicity, denote a± = fii 0A*2> c 2 = /^ 2 0A*i and 

B = g t (BQi, B\ = Qi (Be Qt 
From Proposition 9.1 and the results of Section 8, it follows that 

(9.5) N R2Rl (n) = N ai (n) 
for neN. Now, we need to prove that 

(9.6) N ai (n) = \D 2n (e)\ 

for k e N. Note first that any f-walk w e F(e\) on B\ must begin and terminate with 
an edge of color 1 since no edge of color 2 is incident on e±. Since B\ is obtained from 
£>i by attaching a loop of color 2 to the root e±, the set of rooted alternating d- walks 
on B[ of lenght 2n is in 1-1 correspondence with the set of rooted alternating f-walks 
on Bi of lenght 2n — 1. In fact, each rooted alternating d-walk of lenght 2n on B[ is of 
the form w = (ui,u 2 ), where u\ is a rooted alternating f-walk of lenght 2n — 1 and u 2 
is the rooted loop of color 2. Therefore, 

(9.7) |£>2n(ei)| = |F2n-l(ei)|, 

where D 2n (ei) refers to the loop product and F 2n -i(ei) - to the usual product. This fact 
will be used in the induction proof given below. If n = 1, we use this correspondence 
to get N ai (l) = Ar M1 (l) = |F 1 (e 1 )| = |D 2 (ei)|. Suppose now that (9.6) holds if n is 
replaced byl^/c^n — 1, where o\ is understood to be the s-free convolution of any 
a' and a" (thus, in particular, /i 2 and /ii). Using the proof of Theorem 6.1, we have 

n 

N ai (n) = J] iV m (r) }] N a2 (h)N a2 (k 2 ) . . . N a2 (Av_i). 

r=l fci+fc2 + ...+fc r _i=n— 1 

In view of (9.7), it is enough to show that the expression on the RHS is equal to 
|^2n-i(ei)|. By the inductive assumption, we replace in the above formula each N cr2 (ki) 
by |I?2fei(ei)|, and then, using (9.7), by li^-i^i)!, 1 ^ i ^ r — 1. Therefore, we need to 
justify the above formula, viewing N ai (n) as the number of rooted alternating f-walks 
on B\ and N a2 (fcj), l^i^r — 1, as numbers of rooted alternating f-walks on B 2 . Since 

B 1 =g 1 h B 2 

we can observe that each rooted alternating f-walk of lenght In — 1 on B\ consists of 
a rooted f-walk c = (v ,vi, . . . , v r ) e F r {e{) of color 1, for some 1 ^ r ^ n, with r — 1 
alternating f-walks Wi e F{v,j) on copies of B 2 attached to vertices i>j, where 1 ^ i ^ r— 1. 
The latter must begin and terminate with edges of color 2 since no edge of color 1 is 
incident on the roots of these copies. Thus, the expression given on the RHS of the 
above formula gives the numbers of all alternating rooted f-walks of total lenght 

r + (2*4 - 1) + (2& 2 - 1) + . . . + (2£v_! - 1) = r + 2(n - 1) - (r - 1) = 2n - 1 

summed over 1 ^ r ^ n, which completes the proof. ■ 

Example 9.2. The lowest order first return moments A^ 10M2 (n) are given by 

7W 2 (1) = W w (l)> 
N^ 2 {2) = A^(2)AU1), 
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iW 2 (3) = ^(3)^(1) +^(2)^(1)^(2), 

= ^ w (4)^(1) +2^(3)^(1)^(2)^(1) 
+JV W (2) (i) jv m ( 3 ) + A£ (2) JV m (2) JV w (1) 

where we used the results of Example 6.2. The enumeration of d-walks given in Example 
9.1 can be easily verified by substituting to the above formulas only the non- vanishing 
moments: N^(l) = 1, A^(2) = 2, N^{2) = 2. We have N ai (l) = 1 = D 2 {e x ), 
N n (2) = = D 4 (ei), ^(3) = 4 = £> 6 (ei) and iV CT1 (4) = = D 8 ( ei ), where <n = 

Finally, we will prove the Multiplication Theorem for free independence. In this 
case, one does not need to introduce a new type of graph product since no 'unitization' 
of adjacency matrices is necessary to make them freely independent (basically, this is 
because units are identified in the case of the free product of algebras). Nevertheless, 
if neither of the graphs, Q\ or Q 2 , have loops at their roots, then the set of rooted 
alternating d-walks on Q 1 * Q 2 is empty. 

Theorem 9.2. The Multiplication Theorem holds for free independence, the associated 
product graph Q 1 * Q 2 , and the multiplicative convolution /i X \x\/j> 2 . 

Proof. First of all, we know [2] that 

A(g 1 * g 2 ) = s 1 + s 2 

where S L = 2a°=i A(&), l e I. Moreover, the pairs (Si,S 2 ) and (Si — 1, S 2 — 1) are 
free w.r.t. ip. Also, it is clear that the moments of S 2 Si agree with the moments of 
A*iSa*2- It remains to be shown that the latter coincide with the corresponding numbers 
of rooted alternating d-walks. Using Proposition 2.1, we have \i\ \x\ fj, 2 = ^ Q (/i 2 0/ii) 
and therefore, we can employ the combinatorial formula used in the proof of Theorem 
3.1 to get 

n 

N ai (n) = X>m» 2 N a2 (k 1 )N a2 (k 2 )...N a2 (k r ). 

r=l fci+fc2 + ...+fc r =n 

Now, we can decompose the free product of graphs as 

with the root e of Q\ * Q 2 obtained by identifying the root e± of Q\ with the root e 2 of 
B 2 (see [2,16,20]). In view of Proposition 9.1, 

N C2 (h) = |F 2fci _ 1 (e 2 )|, 

where e 2 is the root in B 2 . Of course, iV Ml (r) = |F r (ei)|, where t\ is treated as the 
root of Qi. Therefore, the RHS of the above combinatorial formula gives the number 
of rooted walks w on Q\ * Q 2 consisting of an f-walk c = (v , v±, . . . , v r ) e F r (e), inter- 
laced with rooted alternating f- walks W\, w 2 , . . . , w r on copies of B 2 attached to vertices 
Vi, v 2 , . . . , v r , respectively. In other words, we have 

w = (p 1 ,w 1 ,(3 2 ,w 2 , . . .,p r ,w r ), 
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Figure 5. An example of the free product 

and therefore, it consists of two rooted alternating f-walks, namely 

ui = (Pi,w 1 ,(3 2 ,W2, ■ ■ ■ ,P r ) and u 2 = w r , 

where Ui begins and terminates with an edge of color 1, and u 2 begins and terminates 
with an edge of color 2. Therefore, w is a rooted alternating d-walk of lenght 2n which 
is 'subordinate' to an f-walk c of lenght r. The summation over 1 ^ r ^ n gives all 
rooted d-walks of lenght 2n, which completes the proof. ■ 

In view of this result, the free product of graphs is 'complete' in our 'category' of 
product graphs since it is naturally related to the rooted alternating d-walks and no 
additional loops are necessary to fit it into the scheme of the Multiplication Theorem. 

Example 9.3. The first return moments iV^p^n) = ^^^(n) of lowest orders 
can be computed using (2.1) and the results of Example 3.2, with the iV^n), where 
o"i = /^2) taken from Example 9.2. For the free product in Fig. 5, we obtain 
D 2 (e) = 0, D 4 (e) = 2, D 6 (e) = and D 8 (e) = 16. 
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